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Abstract: 

We present a systematic study of twist three light-cone distribution amplitudes of vector 
mesons in QCD, which is based on conformal expansion and takes into account meson and 
quark mass corrections. A complete set of distribution amplitudes is constructed for p, 
u, K* and (ft mesons, which satisfies all (exact) equations of motion and constraints from 
conformal expansion. Nonperturbative input parameters are calculated from QCD sum 
rules, including an update of SU(3) breaking corrections in the leading twist distributions. 
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1 Introduction 



The success of QCD as fundamental theory of strong interactions is intimately tied to its 
ability to describe hard inclusive reactions, which has been tested in numerous experiments. 
Progress in the theory and phenomenology of hard exclusive processes has been more 
modest for several reasons. First of all, exclusive processes are in general more difficult to 
study experimentally as they constitute only a small fraction of the inclusive rates at com- 
parable momentum transfers. In addition, there is growing understanding that, although 
the quark counting rules appear to start working at small momentum transfers, QCD fac- 
torization in its standard form |2| may be valid only at very large momentum transfers. This 
is in contrast to inclusive processes like deep-inelastic scattering, where the leading twist 
factorization approximation is adequate already at Q ~ 1 GeV. Evidently, the situation 
calls for a systematic study of preasymtotic corrections to hard exclusive amplitudes. 

Not much is known yet about these preasymptotic corrections. The task is complicated 
by the fact that it actually comprises two different problems. First, it is not known how to 
combine higher twist contributions to hadron form factors with the so-called "soft" or "end- 
point" contributions which are of the same order of magnitude. Second, in order to be able 
to calculate higher twist corrections, one needs to know both the dependence of the leading 
Fock state wave function - with a minimal number of (valence) partons - on transverse 
momentum and the distribution amplitudes with a non-minimal parton configuration with 
additional gluons and/or quark- ant iquark pairs. These two effects are physically different, 
but related to each other by the (exact) QCD equations of motion; taking into account 
one effect and neglecting the other is inconsistent with QCD, unless kinematic suppression 
or enhancement of a particular mechanism can be established. One has to find a basis of 
independent distributions and to develop a meaningful approximation to describe them by 
a minimal number of parameters. 

In this paper we address the second problem only, leaving aside the questions of "end- 
point" contributions and how to generalize the factorization formalism beyond leading twist. 
The main goal of our study is to find out whether higher twist components in hadrons 
have an economic description in QCD. Light-cone distributions beyond leading twist were 
previously addressed in Refs. || £|, [|, || [7|, ||. The existing results are, however, far from 
being complete and sometimes even contradictory. The aim of this paper is to develop 
a systematic formalism for constructing a basis of independent higher twist distributions 
in such a form that it is possible to model them while automatically including all QCD 
constraints. Our approach is an extension of earlier work with the basic idea that the 
equations of motion can be solved order by order in the expansion in conformal spin. These 
relations are exact in perturbation theory although conformal symmetry is broken beyond 
one loop. Taking into account a few low-order terms in the conformal expansion, one obtains 
a consistent set of distribution amplitudes which involve a minimal number of independent 
nonperturbative parameters that can be estimated from QCD sum rules (or, eventually, 
calculated on the lattice). Although in principle this can be done for arbitrary twist, we 
concentrate on twist three distribution amplitudes of vector mesons in this work. A detailed 
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treatment of twist four distributions will be presented elsewhere. 

In Ref . , this program was realized for the pion for which a complete set of distributions 
of twist three and twist four is now available. Vector mesons bring in the complication of 
polarization and are more difficult to treat as the meson mass cannot be neglected. This 
requires a generalization of the techniques of JjJ, which we work out in the present paper. In 
addition, we take into account SU(3) flavour violation effects induced by quark masses and 
construct a complete set of twist three distribution amplitudes for p, uj, K* and <fi mesons, 
which is our main result. All necessary nonperturbative constants are calculated from QCD 
sum rules and the scale-dependence is worked out in leading logarithmic approximation. 

Apart from providing the leading corrections to hard exclusive amplitudes, twist three 
distributions are of special interest as they are free from renormalon ambiguities (power 
divergences of the corresponding operators, in a different language) and their evolution 
with Q 2 is simple in the limit of a large number of colours, N c — > oo, as will be clarified 
in this work. One may hope that these distributions are accessible experimentally. Some 
immediate applications of our results, which we do not pursue in this paper, are to exclusive 
semileptonic and radiative B decays and to hard electroproduction of vector mesons at 
HERA. 

The paper is organized as follows. Section 2 is mainly introductory. We collect neces- 
sary definitions and explain basic ideas. Section 3 contains a detailed study on chiral-odd 
distribution amplitudes, including solution of the equations of motion, conformal expansion 
and renormalization. A similar program is carried out for chiral-even distributions in Sec. 4. 
Section 5 contains explicit models for the p, K* and <fi meson distribution amplitudes up 
to twist three, which involve a minimal number of nonperturbative parameters and satisfy 
all QCD constraints. The final Sec. 6 is reserved for a summary and conclusions. The 
paper contains three appendices: in App. A we collect some useful formulae about orthog- 
onal polynomials, in App. B we elaborate on the structure of conformal expansion for the 
so-called Wandzura-Wilczek contributions, and App. C contains QCD sum rules for the 
nonperturbative expansion coefficients used in Sec. 5. 

2 General Framework 
2.1 Kinematics and notations 

Amplitudes of light-cone dominated processes involving vector mesons can be expressed 
in terms of matrix elements of gauge-invariant nonlocal operators sandwiched between the 
vacuum and the vector meson state, 

(0\u(x)T[x,-x}d(-x)\p-(P,X)}, (2.1) 

where T is a generic Dirac matrix structure and where we use the notation [x, y] for the 
path-ordered gauge factor along the straight line connecting the points x and y: 

[x, y] = Pexp[igfdt (x - y)^{tx + (1 - t)y)\. (2.2) 
J o 
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To simplify notations, we will explicitly consider charged p mesons; the distribution ampli- 
tudes of p° and of K* and mesons can be obtained by choosing appropriate SU(3) currents. 
In order to be able to study SU(3) breaking effects, we keep all quark mass terms. 

The asymptotic expansion of exclusive amplitudes in powers of large momentum transfer 
is governed by contributions from small transverse separations between the constituents, 
which are obtained by expanding amplitudes like ( |2.1| ) in powers of the deviation from 
the light-cone x 2 = 0. To implement the light-cone expansion in a systematic way, it 
is convenient to use light-like vectors. Let P M be the p meson momentum and its 
polarization vector, so that 



P 2 = m% e W. e W = -l, 
We introduce light-like vectors p and z with 

p 2 = o, z 2 = 
such that p — > P in the limit m 2 — > and z — > x for x 

1 



P ■ e (A) = 0. 



(2.3) 



0. From this it follows that 



P 



fj, 2 

nip 



xP — J (xP) 2 — x 2 m 2 



x„ - \p„— + 0(x 4 ) 
zp 



2 -n- 



Useful scalar products are 



x, 
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x m p 
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mz 
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2 M 



pz 
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zp 
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(2.4) 



zP 



p ■ e 



zp 

(A) 



xP) 2 — x 2 m 2 



pi 



m 



2pz 



p z-eM. 



(2.5) 



The polarizaton vector can be decomposed into projections onto the two light-like 
vectors and the orthogonal plane: 



3 (A) 



pz 

;W ■ z 
pz 



z) (eW-p) (A) 



2 



7?) 



2p^ 



;2.6) 
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Note that 

( e W . z ) = ( e M -x). (2.7) 
In terms of the original variables one has 



PfiixP) - x^m 



2 



JX) - ( r M T y^ J ' J - l Wp I r W coal 

We also need the projector onto the directions orthogonal to p and z, 

giu = 9nv - — (p„z„ + p u z^), (2.9) 
pz 

and will often use the notations 

a. = a^, a* = /(pz), (2.10) 

for an arbitrary Lorentz vector a^. 

We use the standard Bjorken-Drell convention 0] for the metric and the Dirac matrices; 
in particular 75 = ij ^ 1 ^ 2 ^ 3 , and the Levi-Civita tensor e^ v \ a is defined as the totally 
antisymmetric tensor with €0123 = 1- The covariant derivative is defined as = D ^ = 
d^ — igA^, which is consistent with the gauge phase factor Q2.2|) , and we also use the notation 
Dfj, — 9 n + igA^x) in later sections. The dual gluon field strength tensor is defined as 
Gnu = ^ fJl upaG pa . Finally, we use a covariant normalization for one-particle states, i.e. 
(p-(P,\)\p-(P',\>)) = (2n) 3 2P°5W(P-P')5 xx , 



2.2 Classification of two-particle distribution amplitudes 

By definition, light-cone distribution amplitudes are obtained from Bethe-Salpeter wave 
functions at (almost) zero transverse separations of the constituents, 

± d 2 k ± 4>(x,k ± ) , (2.11) 

and are given by vacuum-to-meson matrix elements of nonlocal operators on the light-cone: 

(0\u(z)V[z,-z]d(-z)\p-(P,X)). (2.12) 



Note that unlike in Eq. (|2.1|) the separation between the quark and the antiquark is strictly 
light-like. The expansion of ( |2.1| ) near the light-cone in terms of operators with light-like 
separations is the subject of the operator product expansion and considered at length in 
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It turns out that all quark- antiquark distribution amplitudes of vector mesons can be 
classified in the same way as the more familiar nucleon structure functions (parton distri- 
butions) which correspond to the independent tensor structures in matrix elements of type 
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Twist 


2 3 4 

0(1) o(i/Q) o(i/g 2 ) 


spin ave. 

Si 

s± 


fx e U 
9i h± 9s 
hi gr h 



Table 1: Spin, twist and chiral classification of the nucleon structure functions. 



(N(P, S)\if)(z)r[z, —z]ip(—z)\N(P, S)) over nucleon states \N(P, S)) with momentum P and 
spin S. 



Jaffe and Ji |TT] find nine independent quark distributions whose spin, twist and chiral 
classifications are shown in Tab. |l|. The parton distributions in the first row are spin- 
independent, those in the second and third rows describe longitudinally (S\\) and transversely 
(S±) polarized nucleons. Each column refers to twist: a distribution of twist t contributes 
to inclusive cross sections with coefficients which contain t — 2 or more powers of l/Q. The 
underlined distributions are referred to as chiral-odd, because they correspond to chirality- 
violating Dirac matrix structures T = {cx^rys, 1}. The other distributions are termed chiral- 
even, because they are obtained for the chirality- conserving structures T = {7^, 7^75}. 

The nucleon parton distributions are defined by 

(N(P,SMz)a^ l5 [z,-z]^-z)\N(P,S)) 



(S±»p v - S ±u p ll ) I dxe 2ixp - z hx(x,fi 2 



+ {Si_uZ v — S± u z ; 



S ■ z 

^Jp^zf 

M 2 



M 



p ■ z J -1 



dxe 2ixp - z h L {x,fi 2 ) 
dxe 2ixp - z h 3 (x,fx 2 ) 



(N(P,S)\ip(z)[z,-z)ip(-z)\N(P,S)) =2M J dxe 2ixp - z e(x, /x 2 ), 



(N(P,S)\^z)^[z,-z]^-z)\N(P,S)) 
= 2 

(N(P,S)\i>(z) W [z,-zfy{-z)\N(P,S)) 



M 2 r 1 

pj dxe 2lxp - z h{x^ 2 ) + Zll / dxe 2ixp - z f 4 (x,ti 2 ) 

' 1 p ■ z J-l 



(2.13) 
(2.14) 



(2.15) 
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o(i/g) 


o(i/g 2 ) 


e ll 


h 




#3 


e_L 


0± 


#1 , #1 





Table 2: Spin, twist and chiral classification of the p meson distribution amplitudes. 



2M 



S-Z fl 



+ ^ 



5-2 

• 



M 2 J 1 dxe 2ixp - Z g 3 (x,fx 2 ) 



(2.16) 



where P M and 5^ are decomposed similarly to ( |2.4| ) and (|2.6|), respectively, with m p replaced 
by the nucleon mass M; x is the Bjorken scaling variable. The nucleon spin vector is 
normalized as S ■ S = — 1, which differs from the definition used in by a factor M 2 . 

The analysis of vector meson distribution amplitudes reveals an analogous pattern, which 
is no surprise as the operator structures are the same and the p meson polarization vec- 
tor formally substitutes the nucleon spin vector in the Lorentz structures. We find eight 
independent two-particle distributions whose classification with respect to spin, twist and 
chirality is summarized in Tab. 0. One distribution amplitude is obtained for longitudinally 
(en) and transversely (ej_) polarized p mesons of twist 2 and twist 4, respectively. On the 
other hand, the number of twist 3 distribution amplitudes is doubled for each polarization. 
Due to this analogous structure we take over the notations from Tab. |l]for some quantities. 
Again, the higher twist distribution amplitudes contribute to a hard exclusive amplitude 
with additional powers of 1/Q compared to the leading twist 2 ones. The underlined dis- 
tribution amplitudes are chiral-odd, the others chiral-even. Because the matrix element 
(0\u(z)T[z,—z]d(—z)\p~(P,X)) depends on the polarization vector linearly, there is no 
spin-independent distribution amplitude. This is in contrast to the nucleon parton distri- 
butions of Tab. [I], where the dependence on the polarization vector is determined by the 
density matrix (1 + r y^)/2 which contains a spin-independent part. One more difference is 
that e^> for the p meson is a vector, while for the nucleon is a pseudovector. Thus, an 
insertion of additional 275 is necessary in order that matrix elements of relevant operators 
have the analogous Lorentz decomposition. 

The explicit definitions of the chiral-odd p distributions are: 

(0\u(z)a fJ , I/ [z, -z]d(—z)\p~ (P, X)} — 



[efyh,-e$pj / due^^p 2 ) 
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(p-z) 2 p Jo 



(0|u(z)[*, -z]d(-«)|p"(P, A)) = -i (fj - f p 



p ■ z Jo 

t _ j m u + m d 
m„ 



due^ p - z h 3 (u,ii 2 ) 



(2.17) 



p Jo 11 

(2.181 



while the chiral-even distributions are defined as (cf. 
(0\u(z) lp [z,-z]d(-z)\p-(P,\)) = 



f P m p 



p, / due^^^+e^J due^ z g { l\u,p 2 ) 

p ■ z Jo Jo 



1 ■ z r 1 
-z„- -m 2 / due^ pz g 3 (u,p 2 ) 



2 M (p-z) 2 



o 



(2.19) 



and 



(0\u(z) W [z,-z]d(-z)\p~(P,X)) 



o I fp fp 



,m u + m d 



m p e;^ef v p a z p [ due^ z gf '(u, /x 2 ). (2.20) 



Here and below we use the shorthand notation 

£ = u - (1 - u ) = 2u- 1. 

The distribution amplitudes are dimensionless functions of u and describe the probability 
amplitudes to find the p in a state with minimal number of constituents — quark and 
antiquark — which carry momentum fractions u (quark) and 1 — u (antiquark), respectively, 
and have a small transverse separation of order The nonlocal operators on the left- 

hand side are renormalized at scale /i, so that the distribution amplitudes depend on p as 
well. This dependence can be calculated in perturbative QCD and will be considered below 
in Sees. 3 and 4. 

The vector and tensor decay constants f p and fj are defined as usually as 

(0|«(0) 7M d(0)|p-(P,A)) = f p m p ef, (2.21) 
(0|«(0)^d(0)|p-(P,A)) = itfie^Pv-eMp^ (2.22) 



and their numerical values are discussed in Sec. 5. 



7 



All eight distributions (f) = {<f)\\, <p±, g±\g± \ hffl , hi , hs, g%} are normalized as 







~du<f)(u) = l, (2.23) 



which can be checked by comparing both sides of the defining equations in the limit z^ — > 
and using the equations of motion. The rationale for keeping the (tiny) corrections propor- 
tional to the u and d quark masses m u and is that it will allow us to calculate the SU(3) 
breaking corrections for K* and mesons. 

Note that the meson-to- vacuum matrix element vanishes for V = 275, because it is not 
possible to construct a pseudoscalar quantity from z^, and e£\ On the other hand, 
(p.20|) would correspond to (|2.15| ) which defines the spin-averaged nucleon distributions and 
is the only exception to the complete analogy between the nucleon distribution functions 
and the p meson distribution amplitudes. In this case it is the difference in parity between 
and Sfj,, which leads to a completely different decomposition of the matrix elements. 

Because of the analogous structure, the twist classification of the various distributions 



does not require a separate study and can be inferred directly from JTTJ] , see Tab. 0. Its 
physical interpretation, however, deserves a discussion. 

One convenient way to understand the twist classification of distribution amplitudes 
directly from their definitions is to go over to the infinite momentum frame. For definiteness, 
we assume for the moment that p is moving in the positive e 3 direction and p + and z~ are 
the only nonzero component of p and z, respectively. Then the infinite momentum frame 



can be understood as p + ~ Q — > 00 with a fixed p ■ z ~ 1 in (|2.17|) - (|2.20|) . From (|2.6Q it 



follows that in this frame (e^ • z) ~ 1 and ~ 1. This determines the power counting in 
Q of all terms in the right-hand sides of ( |2.17| )-( |2.20| ). The first, second, and third terms in 
( [2.17|) and Q2.19Q behave as 0(Q), 0(1) and 0(1/ Q), respectively, and the right hand sides 
of flTI8D and are of 0(1). 

A mathematically similar, but conceptually different approach to twist counting is based 
on the light-cone quantization formalism [[T3], [I], |TT|. In this approach quark fields are 
decomposed into "good" and "bad" components, so that ip = ip++4>- with ip + = (1/2)7*7. i> 



and ^_ = (1/2)7. 7*0- ^ s discussed in ITTJ, a "bad" component ip- introduces one unit of 
twist. Therefore, a quark-antiquark operator of type ud contains twist 2 (u + d + ), twist 3 
(u + d-, u_d + ), and twist 4 (u-dJ) contributions. This explains why the number of twist 3 
distribution amplitudes is doubled compared with the twist 2 and twist 4 ones (see Tab. 0). 

The physical content of this classification is that a "good" component ip + represents an 
independent degree of freedom corresponding to the particle content of the "Fock state". 
On the other hand, the "bad" components are not dynamically independent, but can be 
expressed in terms of the higher components in the Fock wave function with a larger number 
of constituents, in particular corresponding to a coherent quark-gluon pair. Only the twist 2 
distribution amplitudes correspond to the valence quark-antiquark component in the p me- 
son wave function, while the higher twist amplitudes involve contributions of multi-particle 
states. This point will be discussed in detail in Sees. 3 and 4. 

One important comment is in order. The definition of twist based on power counting in 
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the infinite momentum frame is convenient, because it is directly related to the power ofl/Q, 
with which the corresponding distributions appear in the physical scattering amplitudes, 



and hence is frequently employed in recent works ||11|| . On the other hand, this definition is 
not Lorentz invariant and does not match the conventional and more consistent definition 
of twist as "dimension minus spin" of the relevant operators. For example, the nucleon 
structure function which is identified as twist 3 by power counting in fact contains 
contributions of both operators of twist 2 and twist 3. Similarly, the distribution amplitudes 
of vector mesons that were identified as twist 3 above, actually contain contributions of 
twist 2 operators as well. In Sees. 3 and 4 we will study in detail the operator structure of 
twist 3 distribution amplitudes based on the operator product expansion. In this context 
the conventional "operator" definition of twist will be more adequate. The mismatch of 
different definitions of twist has to be kept in mind, but hopefully will not yield confusion. 
To summarize, Eqs. ( |2.17 )-( p.20 ) define a complete set of valence light-cone distribution 



amplitudes and provide full information on the quark- ant iquark component of the Fock 
wave function of the p meson at zero transverse separation. As mentioned above, not 
all of these distributions are independent. In the following sections we will derive exact 
relations between the twist 3 quark- ant iquark distribution amplitudes and those involving 
one additional gluon, which are introduced below. 

2.3 Three-particle distribution amplitudes of twist three 

Higher Fock components of the meson wave function are described by multi-particle distri- 
bution amplitudes. In this paper we will explicitly deal with three-particle twist 3 quark- 
ant iquark-gluon distributions, defined as 

(0\u(z)-f a [z } vz]gG^(vz)[vz, —z]d(—z)\p~(P, A)) = ip a fr M e^ -p w e^/JV(u,pz) + . . . 

(2.24) 

(0\u(z)-f a j 5 [z,vz]gG^(vz)[vz : -z]d(-z)\p~(P, A)) = p a \pue^-p^&S,]f^A(v,pz) + . . . 

(2.25) 

(0\u(z)a af3 [z,vz]gG^(vz)[vz, -z)d(—z)\p~(P,X)) = 
e (A) . z 

= 2 (p . z) teaPrfP" ~ PPP»9av ~ VaPvQpn + PpPug a n\fIp™>pT{v, pz) + . . . , (2.26) 

where the ellipses stand for Lorentz structures of twist higher than three and where we 
used the following shorthand notation for the integrals defining three-particle distribution 
amplitudes: 

F(v,pz) = J Vae-* pz ^- a « +va ^F(a d ,a u ,a g ). (2.27) 

Here T = {V, A, T} refers in an obvious way to the vector, axial- vector and tensor distribu- 
tions, a is the set of three momentum fractions: a d (d quark), a u (u quark) and a g (gluon), 
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and the integration measure is defined as 

J Va = J da d J da u J da g 5(1 — c^). (2.28) 
The normalization constants /Z, , f£ p are defined in such a way that 

J (a d - a u ) V{a d , a u , a g ) = 1, 

Jl)aA(a d ,a u ,a g ) = 1, 

J Va (a d - a u ) T(a d} a Ul a g ) = 1. (2.29) 

Choosing the normalization in this way we anticipate that the function A is symmetric and 
the functions V and T are antisymmetric under the interchange a u <-> a d in the SU(3) 
limit (cf. 0), which follows from the behaviour of the corresponding matrix elements under 
G-parity transformations. Numerical values of all constants are discussed in Sec. 5. 

With these major definitions at hand, we now proceed to a systematic study of the 
twist 3 distribution amplitudes. 



3 Chiral-odd Distribution Amplitudes 

This section is devoted to the general discussion of chiral-odd distributions of twist 3. We 
demonstrate that the two-particle distribution amplitudes h[?\u, jj 2 ) and h^\u, jj, 2 ) can be 
eliminated in favour of independent dynamical degrees of freedom and expressed in terms of 
leading twist two- and three-particle distributions. The corresponding relations are worked 
out in detail and solved explicitly. A similar relation between the nucleon structure function 



gi and the twist 2 part of is known as Wandzura-Wilczek relation [13]. We also investigate 
the expansion of twist 3 distributions in terms of matrix elements of conformal operators. 
We demonstrate that the equations of motion are satisfied order by order in the conformal 
expansion which provides, for this reason, a systematic approach to the construction of 
models of distribution amplitudes, consistent with QCD constraints. The renormalization 
of all distributions is worked out in the leading logarithmic approximation. 



3.1 Equations of motion 

The basis of twist 3 distributions defined in Sec. 2 is overcomplete. Due to the QCD 
equations of motion, the number of independent degrees of freedom is less than the number 
of independent Lorentz structures, and our first task will be to reveal the corresponding 
constraints. 

The standard technique for this purpose is to derive relations between towers of local 
operators which arise in the Taylor expansion of the nonlocal operators in Eqs. ( |2.17| ), ( [2.18| ) 
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and whose matrix elements are just moments of the distribution amplitudes. A more elegant 
and economic approach is to use exact operator identities between the nonlocal operators 
H [7j (see also In the present context, we need the identities 



_d_ 

dx. 



{u(x)a^ u x v [x, -x]d(-x)} 



= i J dvv u(x)x a a al3 [x,vx}x^gG fl j3(vx)[vx } —x]d(—x) 

— ix^dp {u(x) [x, — x}d(— x)} — (m u — md)u(x) $ [x, — x\d{— x) 
u(x)[x, — x]d(— x) — u(0)d(0) = 



1 ,t 

dt I dv u(tx)x a a a /3[tx,vx]x^gGn l 3(vx)[vx, —tx\d(—tx) 
o J-t 



+ i f dtd a [u(tx)a a/ 3X l3 [tx, —tx}d(—tx)} 



(3.1) 



+i(m u + md) / dt u(tx) $ [tx, — tx]d(— tx). (3.2) 
Jo 

Here we introduced a shorthand notation for the derivative over the total translation: 



d 

d a {u(tx)T[tx, — tx\d(— tx)} = -q-^ {u(tx + y)T[tx + y, —tx + y]d(—tx + y)} 



, (3.3) 



y^O 



with the generic Dirac matrix structure T. 

In the light-cone limit x 2 — > matrix elements of the operators on both sides of Eqs. ( p.l|) 
and Q3.2| ), sandwiched between the vacuum and the p meson state, can be expressed in terms 
of the distribution amplitudes defined in Sec. 2: 

(0\u(x)a f , l/ x u [x, -x]d(-x)\p~ (P, A)) = 
' ' fe (A) • x x 



if, 



(A) 



_ m A_f^l( Xfi _^_ p \ Jj ue i^ h f^^)^^^)A } (3.4) 



Px 



\u(x)[x, —x]d(—x)\p (P, A)) = 

m u + m d \ (e (A) .^2 f 1 du SP--\ h f^^) +0 (x 2 % (3.5) 



-i[fj-fp 
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and likewise for three-particle distributions. 

The matrix elements of Eqs. ([3J]) and ( |3.2|) yield a system of integral equations between 
two- and three-particle light-cone distribution amplitudes:[] 

-ipz £due^ pz ^h^\u) -2 j\ue* pz {h ( ^{u) - <f> ± (u)) = 

= Clp&Y f dvi}T(v,pz) + (1 - 5+) (pz?f Q e* pz h\\u) + i5.pz j\u e^ pz <P\\{u) , (3.6) 
(1 - 5 + ) [ due^ pz h\f\u) = iCLpz [\dt f 1 dvT{v,tpz) + Cdt f due mpz hf{u) 



11 p Jo J-l Jo Jo 

-5+ J^dt £due^ tpz ^\(u), (3.7) 
where we have discarded all corrections of order x 2 , set x^ = z^ and introduced the notations 



, f P m u ±m d T g p 

fp m P fpm p 



Eqs. ( |3"1)| ) and ( |3.7| ) are exact in QCD. Note the terms with total derivatives in ( |3.1| ) and 
(P-2|), which induce mixing between h^\u) and h f 1 (u). Such contributions are specific for 
exclusive processes and have no analogue in deep-inelastic scattering. Note also that quark 
mass corrections bring in the leading twist chiral-even distribution <p\\(u). 

We can solve Eqs. ( |3.(j| ) and ( ^.7[ ) for (u) and h!fi (u) in terms of the other distributions. 
To simplify the algebra, it is convenient to consider moments in an intermediate step. 
Defining 

M^ = £duCHA^l M^= £duCh^'®(u) (3.9) 



and 



T=0 



J (a d -a u - va g ) n T(a d , a u , a g ) (3.10) 
and expanding ( |3.6|) and Q3.7|) in powers of (pz), we obtain 

(3.11) 

(1 - 6+) M« -^—M^+5 + ^— M\ - f dv T n ^(v) = 0. (3.12) 

7i+l Ti + 1 'n + lJ-i 



lr rhe suppressed corrections 0{x 2 ) drop out; the 0(x 2 ) term in the Lorentz structure in the second line 



of (3.4), however, does give a contribution to the left-hand side of (3.1) after taking the derivative with 
respect to x M . 
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Combining these two equations, one gets the following recurrence relations for and h\^\ 
(n + 2) M® - n M® 2 = 2M^ + (J p f_dv {(n - l)n vT n _ 2 (v) + (n - 2)n T n _ 3 (v)} 

-S+nM^-S-nM^, (3.13) 

(1 - 6+) {(n + l)(n + 2)M« - (n - \)nM$_ 2 } = 

= 2M^ + Clpf^ dv {n{n + 2) T n . x {v) + (n - l)nvT n _ 2 {v)} 

-5 + {n + 2)Ml-5.nMl_ 1 . (3.14) 

Recurrence relations of this type are easily solved by transforming them into differential 
equations. For instance, for the distribution amplitude h^\u) one finds a second order 
equation: 

(l-6 + )uu(h\? ) )''(u) = -$(u) (3.15) 



with 



(«) = 2Mu) ~ S + (0|| H - + \8-<t>\\{u) 

+ ( Sp~r [da d [da u - U/- + a u -^ 1 J T(a). (3.16) 

era J J 1 — a u — ad \ dad da u ) 



Here and below we use the shorthand notation u — 1 — u. The solution of this equation 
with boundary conditions specified by the values of the first two moments reads 

r 1 } 1 

(1 - 5+) h\°\u) — u dv- $(u) + u dv- . (3.17) 

' « 

The solution for (u) can be obtained in a similar manner and reads: 
h®(u) = \t(jdv±*(v) - fdv±*(v)\ +M||(«) 



u it 



+ Cj P T" [da d [da u - T(a). (3.18) 

h du J J 1 — a u — a d 



o o 
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According to the various "source" terms on the right-hand side of ( |3.13| ) and fl3.14| ), one 
can decompose the solution in an obvious way into three pieces as 



hf{u) = hf WW {u) + hf g {u) + hf m {u), (3.19) 

hf{u) = h\( )WW (u) + h\° )9 (u) + h\[ )m {u), (3.20) 

where n» W (u) and h^ WW (u) denote the "Wandzura-Wilczek" type contributions of 
twist 2 operators, h^ g (u) and h^ 9 (u) stand for contributions of three-particle distribu- 
tions and hfi m (u) and h^' m (u) are due to the quark mass corrections. In particular, we 
get 

hf WW {u) = t ( Pdv^l _ [ l dv <tM\ , (3.21) 



hP WW (u) = 2 (« f dv^ + u I 1 dv*±M) . (3.22) 



o 



These are the analogues of the Wandzura-Wilczek contributions to the nucleon structure 
functions gr{x,Q 2 ) H and h L (x,Q 2 ) [[□]]. 



The relations Eqs. ( |3.17| ) and ( |3.18|) are the principal results of this section: chiral-odd 
two-particle distribution amplitudes of twist 3 are expressed in terms of the leading twist 
amplitudes and the three-particle twist 3 distribution. In the next subsection we will discuss 
how to proceed further with this rather complicated formal solution, concentrating on the 
massless quark limit. 



3.2 Conformal expansion 

The conformal expansion of light-cone distribution amplitudes is analogous to the partial 
wave expansion of wave functions in standard quantum mechanics. In conformal expansion, 
the invariance of massless QCD under conformal transformations substitutes the rotational 
symmetry in quantum mechanics. In quantum mechanics, the purpose of partial wave 
decomposition is to separate angular degrees of freedom from radial ones (for spherically 
symmetric potentials). All dependence on the angular coordinates is included in spherical 
harmonics which form an irreducible representation of the group 0(3), and the dependence 
on the single remaining radial coordinate is governed by a one- dimensional Schrodinger equa- 
tion. Similarly, the conformal expansion of distribution amplitudes in QCD aims to separate 
longitudinal degrees of freedom from transverse ones. All dependence on the longitudinal 
momentum fractions is included in terms of functions (orthogonal polynomials) forming irre- 
ducible representations of the so-called collinear subgroup of the conformal group, SL(2,R), 
describing Mobius transformations on the light-cone. The transverse-momentum depen- 
dence (the scale-dependence) is governed by simple renormalization group equations: the 
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different partial waves, labelled by different "conformal spins", behave independently and 
do not mix with each other. Since the conformal invariance of QCD is broken by quantum 
corrections, mixing of different terms of the conformal expansion is only absent to leading 
logarithmic accuracy. Still, conformal spin is a good quantum number in hard processes, 
up to small corrections of order a 2 s . Application of conformal symmetry to the studies of 
exclusive processes in leading twist have received a lot of attention in the literature, see 
e.g. pi 0, 0. 

Despite certain complications, the conformal expansion presents a natural approach to 
the study of higher twist distributions, which has even more power than in leading twist. 
The reason is that conformal transformations commute with the exact QCD equations of 
motion since the latter are not renormalized^. Thus, the equations of motion can be solved 
order by order in the conformal expansion. In this section, we use the approach of [|l^. [7| 
to work out the explicit form of the conformal expansion for the chiral-odd distributions 
(f)±(u),h^\u),h^\u) and T(a), and solve the constraints (|3.19|) and (|3.20| ) order by order 
in conformal spin. 

Since quark mass terms break the conformal symmetry of the QCD Lagrangian explicitly, 
one might expect difficulties to incorporate SU(3) breaking corrections in the formalism. In 
fact, the inclusion of quark mass corrections turns out to be straightforward and produces 
two types of effects. First, matrix elements of conformal operators are modified and in 
general do not have the symmetry of the massless theory. This is not a "problem" , since the 
conformal expansion is designed to simplify the transverse momentum dependence of the 
wave functions by relating it to the scale dependence of the relevant operators. This depen- 
dence is given by operator anomalous dimensions which are not affected by quark masses, 
provided they are smaller than the scales involved. Second, new higher twist operators arise, 
in which quark masses multiply operators of lower twist, see the previous section. These 
additional operators, again, do not pose a "problem" and can be expanded systematically 
in conformal partial waves, leaving the quark masses as multiplicative factors. Explicit ex- 
amples are considered later in Sec. 5, while in this section we neglect operators proportional 
to the quark masses and set S± = in the formulae obtained in Sec. 3.1. 

The conformal expansion of distribution amplitudes is especially simple when each con- 
stituent field has fixed (Lorentz) spin projection onto the light-cone. Such constituent fields 
correspond to the so-called primary fields in conformal field theories, and their conformal 
spin equals 

j = 1(1 + 8), (3.23) 

where / is the canonical dimension and s the (Lorentz) spin projection. Multi-particle states 
built of primary fields can be expanded in increasing conformal spin: the lowest possible spin 
equals to the sum of spins of constituents, and its "wave function" is given by the product 
of one-particle states. This state is nondegenerate and cannot mix with other states because 
of conformal symmetry. Its evolution is given, therefore, by a simple renormalization group 
equation and one can check (see Sec. 3.3) that the corresponding anomalous dimension is 

2 More precisely, one can regularize the theory in such a way as to preserve the equations of motion. 
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the lowest one in the whole spectrum. Therefore, this state is the only one which survives in 
the formal limit Q 2 — > oo; following established tradition we will refer to the multi-particle 
state with the minimal conformal spin as "asymptotic distribution amplitude" . 

An explicit expression for the asymptotic distribution amplitude of a multi-particle state 
built of primary fields was obtained in Refs. [IS, 0]: 



T[2jiJ . ..T[2j m \ 

Here the jj. are the conformal spins of the constituent fields (quark or gluons with fixed 
spin projections). This distribution has conformal spin j = ji + ■ ■ ■ + j m . Multi-particle 
irreducible representations with higher spin j + n,n = 1,2,... are given by orthogonal 
polynomials of m variables (with the constraint Y^k=i a k — 1 ) with the weight function 
(PI- 

A classical example is the leading twist quark- ant iquark distribution amplitude. The 
distribution amplitude <fix(u), defined in ( |2.17| ), has the expansion 

oo 

^(u) = 6uu£aiCf , (fl, (3.25) 

n=0 

where C%/ 2 (£) are Gegenbauer polynomials (see e.g. jnj). The dimension of quark fields is 
I = 3/2 and the leading twist distribution corresponds to positive spin projection s = +1/2 
for both the quark and the antiquark. Thus, according to ( p.23|) , the conformal spin of each 



field is j q = jg = 1; the asymptotic distribution amplitude ( |3.24| ) equals 4> as (a q} otq) = Qa q aq 
and has conformal spin j = 2. Taking into account a q + otq = 1 and denoting u = a q we 
arrive at the first term in the expansion ( |3.25| ). The Gegenbauer polynomials correspond to 
contributions with higher conformal spin j + n and are orthogonal over the weight function 
Quu. 

Note that = 1 due to the normalization condition ( |2.23| ). In the strict massless 
limit only the terms with even n survive in Eq. (|3.25|) because of G-parity invariance. The 



conformal expansion, however, can be performed at the operator level and is disconnected 
from particular symmetries of states such as G-parity. The following discussion is, therefore, 
valid for arbitrary n. We keep terms with n = 2k+l for the later discussion of SU(3) breaking 
corrections. 

The conformal expansion of the twist 3 two-particle distribution amplitudes (u) and 

h\f\u) is less immediate. As a first step, one has to decompose them into components 
built of primary fields - - with fixed spin projections. To this end, we define a set of 
auxiliary amplitudes h'^(u) and h^(u) using the spin projection operators P + = (1/2)7*7. 
and P_ = (1/2)77* to single out quark states with s = +1/2 and s = —1/2, respectively, 
(see 0,|3): 

(0\u(z)w[z,-z}d(-z)\p-(P,\)) = fJml——j o due^h^iu), (3.26) 
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(0\u(z)^ r [z,-z]d(-z)\p-(P,X)) = f p ml ( —^ f'due^h^iu), (3.27) 

H p ■ z Jo 

which are related to h^\u) and h\°\u) by (see ( PT7| ), (gig)) 



II 

hU( tt ) = /l W( u ) + iZX ; p j (3.28) 

/i (it) = -hf(u) + (3.29) 

11 2 cm 

The conformal expansion of (it) and h^* (it) is straightforward and is given by 

oo 

hK(u)=2*E>>M VS >®, (3-30) 

n=0 

oo 

^( U )=2«5:^P ? (o,i) (0) (3 . 31) 

n=0 

where -Pi°'^(£) are Jacobi polynomials (see e.g. [[HJ) and the n-th term corresponds to 
conformal spin j = n + 3/2. Substituting these expansions in (|3.28|) and (|3.29| ) and using 
the identities ( PO| ), (|A.10|) and ( [7011) in App. we obtain 



fcjjV) = £ (H n -# n _i) Cy 2 (0+ E (^-^-1)^(0, (3.32) 

n=0,2,4,... n=l,3,5,... 

^ = 4,16 [mJ^^^^J^^^ 1 ■ (3 - 33) 

where H_i = h-\ = is implied, and 

- (-!)»/# 



2 



= " + l 2 j " , (3.34) 

for n = 0,1,2,... correspond to G-parity conserving and G-parity violating contributions, 
respectively. Note that the coefficient in front of each orthogonal polynomial in ( |3.32|) and 



33|) does not correspond to a definite conformal spin; in contrast to ( |3.25| ), ( |3.30| ), and 



( |3.31|) , it is rather given by difference between the contributions of two successive conformal 



spins. 

The conformal expansion of the twist 3 three-particle distribution amplitude gives yet 
another example for the general expression Eq. ( |3.24| ). The expansion reads 

oo 

T(atd, a u , 1 - a d - a u ) = 3Q0a d a u (l - a d - a u ) 2 E ^I,iJk,i( a d, ««), (3.35) 

k,l=0 
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where Jk,i(otd, &u) = Jk,i(6, 2, 2, a<i, a u ) are particular Appell polynomials of two variables 
(see p269 of []19|). The conformal spin of a generic term in this expression equals j = 
k + I + 7/2 and is the same for all contributions with equal sum n = k + I. This illustrates 
that three-particle conformal representations are degenerate; the number of independent 
operators with the same spin in fact increases with the spin. Conformal symmetry does not 
allow mixing between contributions with different j = n+7/2; it does allow, however, mixing 
with each other of different states with the same value of j. Therefore, the mixing matrix 
for higher twist operators becomes only block-diagonal in the conformal basis, instead of 
being diagonalized like in leading twist. 

Next, we are going to demonstrate that conformal expansion is fully consistent with the 
equations of motion. To this end we need to show that the conformal expansion coefficients 
for two-particle twist 3 distributions can be expressed in terms of the expansion coefficients 
for three-particle distributions with the same conformal spin, and we need to separate the 
Wandzura-Wilczek contributions. The calculation is straightforward, although somewhat 
tedious. 

We decompose 



H n = H, 



WW 



+ H 9 



h n = h^ w +hl (3.36) 
and start with the Wandzura-Wilczek contributions (|3.21|) and (|3.22|) which give rise to 



auxiliary amplitudes h^ WW (u) with 



ww , 



u 



2u 



dv 
o v 



(3.37) 



h^ WW (u) = 2u 



(3.38) 



The integrals on the right-hand side of ( |3.37| ) and ( |3.38| ) can easily be done using ( |A.3|) and 

(3.39) 



_r*Mo + / i *,M!) = _3£ aiiffl ' tt) . 

n=0 







V 



V 



Substituting the recurrence relations for Jacobi polynomials, (|A.7j) , into this result, one 
immediately obtains 



H 



ww 

n 



H 



WW 



3(w+ 1) 
2n + 3 

3(n + l 



_l. h ww _ 3(ra + 1) | 

a ni ,L n ~ oil U n-1 



2n + 1 



2n + l ' n 2n + 3 ' 



-a 



(n = 0,2,4,...), 



(n = 1,3,5, 



(3.40) 



For even n, we find that which corresponds to the conformal spin n + 2 in the expansion 
for the twist 2 distribution amplitude gives rise to H^ w and which corresponds to the 

conformal spin n + 3/2 and n + 5/2, respectively. Likewise for odd n, gives rise to h^ w 
and hn+[. These values of the conformal spin do not match the expansion in Eqs. ( |3.30[ ) and 
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Q3.31| ). This is, however, not a contradiction since Wandzura-Wilczek terms are in fact not 
intrinsic twist 3 distributions, but correspond to matrix elements of twist 2 operators over 
p mesons with different (longitudinal) polarization. To relate matrix elements of conformal 
operators over longitudinal and transverse p mesons, one has to perform a spin rotation (in 
the p meson rest frame) which does not commute with the generators of collinear conformal 
group. As shown in App. [B], this rotation gives rise to the shift in conformal spin and 
exactly explains the mismatch appearing in Eq. ( |3.40|) . Therefore, the conformal symmetry 
is realized in Wandzura-Wilczek contributions as well, but to see this one has to supplement 
the conformal classification of operators by conformal transformation properties of the meson 
states. 

The three-particle contributions H% and h 9 n of ( [3.36|) can be treated similarly. From 
the solutions for h$ (u) and h*f (u) in Sec. 3.1 we obtain the corresponding auxiliary 
amplitudes: 



u 



u , K(v) 
dv 



where a„ = 1 — ad — a u , and 



K{v) 



+ 



dv 



K{v) 



d z" 1 f u 1 
K(u) = — / da d / da u — 
du Jo Jo a 



+ 



d 

du Jo 



da n 



u i 
da u — T(a) 



du 



/ dad / da u — T(a) 
Jo Jo a„ 



' d d 
a d - h a u - 



da n 



da, 



1 T(a) 



(3.41) 
(3.42) 



(3.43) 



Substitution of the conformal expansion ( |3.35| ) into ( |3.43j ) yields 



^ t k\l\(-l) k 

K(u) = 180 uu 

k,l=0 



k-l 



(1.1) 



(k + l + 2)\ \(k + l + 3) 



(3.44) 



where we have used Eq. QA.15 ) to perform the integration. The final integration involving 
K(v) on the right-hand side of ( |3.41| ) and ( p.42| ) can be done similarly to ( |3.39| ) by using 

: 







v 



v 



180 x: 



k\l\(-l) k+l ( k-l 



CO 



k,l 



(0,0) 



k,l=0 



(fc + Z + 3)! + ^ + 4 



and the last term of (|3.41|) and ( p.42|) can be integrated using (|A.16|) : 



(3.45) 



d_ 

du 



dad / da u — T(a) = 180mm ^ 



k\l\(-l) k 



k-l 



to 



a 



k,l 



Pk+l+2(0 Pk+l+l(0j- 

(3.46) 



(i,i) 



fej=0 



(fc + Z + 3)! \k + l + 3 
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Substituting ( ggg ) and ( ggjj ) into ( PUD , ( ggg ), and using the identities (|Q|)-(|S77|), we 
find 

^ = 180Cs T p E kKn ~t~ 2 o^ + ^ (-ir^iL-^l (- = 3,4,5,...), (3.47) 



fc=0 



n-2 



(n + 2)! 
fc!(n-Jfe-2)!(Jfe + 2) 



^ = -180C 3 T P E-^ 7~~W ' (-l) w "MU-fc-2} (n= 2,3,4,...), (3.48) 

while ifg = iff = iff = hfj = hf = 0. Here we introduced the following quantities: 

T T 

<"jM = (3.49) 

-f^o = ^4^- (3-50) 

We find that the coefficients with fixed k + l = n — 2, which correspond to the conformal 
spin j — n + 3/2, all contribute to H% and h 9 n corresponding to the same conformal spin 
j = n + 3/2, as anticipated. For later convenience, we give the lowest order coefficients 
(ET37D and d^H) explicitly: 



h 9 2 — -15C^wf ,o}' — --^(ip^lifl}! h\ — (J p i - 3^{2,o} )>•••) (3.52) 



where we substitute ^ = 28/3, which follows from the normalization condition ( |2.29| ). 
From (|3.36|) , (|3.40|) , and ( p.47| ) it follows that the two lowest coefficients H and Hi are 



completely determined by the value of cZq = 1, which results in H$ = 1 and H\ = —2. It 
is easy to see that these values for H and Hi ensure Jq 1 du (u) = Jq du (u) = 1 and 
therefore, the normalization condition for <f>±(u) ensures correct normalization of h^\u) and 
h\°\u). 

To summarize, we have demonstrated that the equations of motion that relate different 
twist 3 operators can be solved order by order in the conformal expansion. In other words, 
equations of motion impose "horizontal" relations between operators of the same conformal 
spin and do not involve other spins. This picture is somewhat complicated by the Wandzura- 
Wilczek contributions of the operators of lower (leading) twist which have a more peculiar 
structure. The explicit relations derived above can be made somewhat more compact by 
assuming G-parity invariance. In this case = for odd n, ujJ. 1 = —ujf k , H n = hj^, 
and h n = 0. As mentioned above, the G-parity violating terms are only relevant for SU(3) 
breaking corrections in the distribution amplitudes of K* mesons. 
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3.3 Renormalization and scale-dependence 

The scale-dependence of the chiral-odd distribution amplitudes 0j_(m ;/ u 2 ), h» (u, yU 2 ) and 

h^\u, /i 2 ) is governed by the renormalization group (RG) equation for the relevant nonlocal 
light-cone operators appearing in the definitions (|2.17| ), (|2.18|) , and (|2.26| ). Unlike inclusive 
processes, operators involving total derivatives have to be taken into account since they 
contribute to nonforward matrix elements, and this leads to additional operator mixing. 
For example, consider the set of local operators u(0)(D .) n ~ k a±.(D .) k d(0) {k = 0, 1, 2, . . .), 
which contribute to the n-th moment of the twist 2 distribution amplitude <f)±(u, fi 2 ). These 
operators differ by total derivatives and all mix with each other under renormalization; 
to calculate the scale-dependence, one has to find the eigenvalues and eigenvectors of the 
corresponding anomalous dimension matrix. 

As is well known [15, |TB|, [Ll| , conformal expansion provides the solution to this problem. 



Analysis based on the anomalous Ward identities for the dilatation and special conformal 
transformation (which are members of the conformal group) shows [IS] that, to leading 



logarithmic accuracy, the conformal operators with different conformal spin do not mix 
with each other and thus diagonalize the anomalous dimension matrix. As a result, by 
employing a conformal operator basis we do not encounter any additional operator mixing 
compared to inclusive processes. The relevant anomalous dimensions, which correspond 
to the eigenvalues of the anomalous dimension matrix, can be extracted directly from the 
results for renormalization of the corresponding parton distribution functions. In particular, 
the one-loop anomalous dimensions for hffi(u, jj 2 ) and h^\u,jj 2 ) are the same as for the 
chiral-odd parton distribution functions pj, [22], |23|, |24| as will be shown in the following. 

Our main task in this section is to reveal the explicit operator content of the conformal 
operators, corresponding to particular coefficients in the conformal expansions (|3.32| ), ( |3.33| ), 
( |3.25|) , and ( |3.35| ). We give a one-to-one correspondence between the conformal basis and the 
basis used in the inclusive case. This allows us to determine the anomalous dimensions of the 
conformal operators and to find the evolution of the distribution amplitudes h^\u,fi 2 ) and 

h^\u,fi 2 ) through the conformal expansion. We will work out this program for arbitrary 
conformal spin. 

One complication is that the conformal operator basis for the higher twist operators is 
degenerate (see ( |3.35| )) and the mixing matrix becomes only block-diagonal instead of being 
fully diagonalized like in leading twist. Consequently, the conformal expansion for three- 
particle contributions to hif (u, /x 2 ) and hi (u, fi 2 ) does not resolve possible mixing between 
components with the same conformal spin. This is similar to mixing of the many quark-gluon 
correlation operators for the corresponding twist 3 parton distribution functions [[21], p3f . 
It has been shown recently |25|, ^] , however, that an important simplification occurs in the 
limit of a large number of colors or of large spin (moment of parton distribution function). In 
these limits all complicated mixing disappears and the twist 3 parton distribution functions 
obey simple DGLAP-type evolution. We will demonstrate that the twist 3 distribution 
amplitudes obey a similar pattern. 

Let us start with the Wandzura-Wilczek terms h\^ WW (u, jj, 2 ) and h\\ WW (u, /i 2 ) (see 
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fl3.32| ), ( |3.33| ), and ( |3.36D ). The coefficients H^ w and h^ w in their conformal expansions 



have been expressed by of ( |3.25|) as shown in ( |3.40|) . Therefore, the scale-dependence of 
a^r, i.e. of the twist 2 amplitude 4>±(u,[i 2 ), completely determines that of hfi WW (u, jj, 2 ) and 

h^ WW (u, fi 2 ). From ( [3.25D and orthogonality relations of Gegenbauer polynomials JL9], we 
obtain 

<^ 2 ) = 0/ 2(2 ?,| 3) C duCTmdu^ 2 ). (3.53) 



3(n + l)(n + 2) Jo 
Substituting fl2.17|) into the right-hand side of ( |3.53|) gives 



with 



(^») <"'> = ^l^U)^ ®*^^)) <3 ' 54) 



Sli(x; f ) = (id.)" u(x)ef"o„.Cy ( ^ | d(x), (3.55) 



where the local operator on the right-hand side is renormalized at fi 2 , D— V — D, and 9 M 
is the total derivative ( |3.3j ). fl„ (x, /i 2 ) is the conformal operator of conformal spin j = n + 2 
T6| , |18fPl Therefore, it is RG covariant to leading logarithmic accuracy and satisfies the RG 



equation: 

{^ + m iTg + S 7 ") ^fat*) = °' (3 - 56) 

where 7^" is the one-loop anomalous dimension of the operator Q^. To establish a formal 
connection with the results given in the literature, it is convenient to consider the case 
where fl^ is diagonal in quark flavour corresponding to the flavour matrices A3, Ag, and to 
take the forward matrix element of ( |3.56| ) over the nucleon state \N(P,S)). Because the 
total derivatives drop out in this matrix element, ( p.56|) reduces to the RG equation for 



{N(P } S)\if)(0)a±. {iD.) n x/j(0)\N(P, S)), which gives the n-th moment of the nucleon parton 
distribution function hi(x,/j, 2 ) of (|2.13|) . By matching with the results for renormalization 
of hi(x,fi 2 ) |2l| , we obtain 7,^ as 

7n x = 4CV (V(n + 1) + 1e - ~ + ^) , (3.57) 

where i[)(n + 1) = X^=i — 1e is the digamma function, 7^ is the Euler constant, and 
C F = (N 2 - l)/2N c . From (|3~5l) - (|3~57|) , we obtain 

(fJ^)(Q 2 ) = L^ b (fJa^)(^ (3.58) 



where L = a s (Q 2 )/a s (fi 2 ) and b = (1 1JV C — 2JV»/3. Combined with ( ^32l) , (|Q3D and (|3T40| ) 



this result gives the /^-dependence of h ( fi WW (u, jj 2 ) and h[^ WW (u, jj 2 ) and also determines 



3 In principle, one can construct a tower of conformal operators (d.) m fl^ (to = 0, 1, . . .) with the same 
conformal spin, but with the different "third component" of it. 
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evolution of the twist 2 distribution amplitude (p±(u,fi 2 ) of ( |3.25| ). We note that ( ^.58[ ) for 
n = gives the scale-dependence of the tensor decay constant fj as 



FAQ 2 ) = L c ^f^ 



(3.59) 



because cLq = 1. 

The three-particle contributions h^ 9 (u,fi 2 ) and h^ 9 (u, fi 2 ) can be treated in a similar 
manner, although the discussion becomes more complicated because one has to deal with 
a degenerate representation of the conformal group (see Q3.32j ), ( |3.33| ) and ( |3.36j )). The 
relevant expansion coefficients and h 9 n are expressed by in ( |3.35|) and shown in (ft.47 ) 
and ( |3.48D . Thus the first step is to demonstrate that uj k r l are given by matrix elements of 
the local conformal operators derived in Using (|3.35| ) and the orthogonality relations 
( |A.13| ) for the Appell polynomials, we obtain 



n-2 



Da J k . 



n-k-2 



[oc d , a u )T(a) 



360(-l) n v 
2"+ 1 (n + l)(2n + l)!! ^ 



r,n— r— 2 rv n—r—2,k 



(3.60) 



for k = 0, 1, 



, n — 2, and n — 2, 3, . . . Making use of 

(0| , u(tz)cr' / .[t2;, vz]gG v .(vz)[vz, wz\d{wz)\p~ (P, A)) = 

= (p • z)(e (A) • z)fT p m p J Vae- ip - z(tau+wad+va3) T{a), 
which is equivalent to (|2.26|) , the left-hand side of ( p.60|) gives 

75^ ( eW .,)l P .^- < ' A ^-"' >'"'' p - A »- 



with 



aL(o) 



m k+l j k . 



D v Df 



n(a;)^.^.(0)rf(y) 



(3.61) 



(3.62) 



(3.63) 



x=y=0 



where the covariant derivatives and act on the coordinates x and y, respectively. 
A^ n _ fc _ 2 {k = 0, 1, ... ,n — 2) are the twist 3 conformal operators of spin j = n + 3/2, 
forming a degenerate basis for three-particle representation [18fl.n Inverting the matrix 
W^ l) in ( |3~60D , we obtain 



Nl(~l) k 



90k\(n - k — 2)\ 



where iVj is the dimensionless and scale-independent constant: 



2 n - 1 (n + l)(2n + 1)!! 
m p (e( A ) • z)(p ■ z) 71 ' 1 ' 



(3.64) 



(3.65) 



4 One can generate a tower of conformal operators with a different "third component" of conformal spin 
by acting repeatedly with d. on (3.63). 
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The numerical factor in front of the matrix element ( |3.64j ) is put for later convenience. The 
operators Q kn _ k _ 2 are given by linear combinations of Aj n _ r _ 2 (r = 0, 1, ... ,n — 2) and 
therefore have conformal spin j = n + 3/2. 

The second step is to determine the explicit form of the relevant conformal operator 
®fc,n-fc-2- For the few lowest conformal spins, it is easy to express 0^ n _ fc _ 2 as a linear 
combination of Aj n _ r _ 2 using ( p.60[ )-( p.63[ ), but the procedure becomes more complicated 
for higher conformal spins. Using orthogonality relations for the Appell polynomials, ( |A.14j ), 
it proves possible, however, to determine Q k n _ k _ 2 for general n up to total derivatives. 



/ 



Vaa k d a n u - k - 2 T(a 



360(-l) n £;!(n- jfe-2)! . , 

^kn-k-2 r. T77 7T7^ 7T7T + (terms involving W/__i_2 

' 2" • '(;/ + l)(2n + 1)!! - 

1 



— (0\m(iD.) n - k - 2 o v .9G v .(0)(iD.) k d(0)\p-(P,X)), (3-66) 



/a^ m p( e(A) • z)(p- *) 

where the last line is obtained by substituting ([3.61 ) into the left-hand side. From (|3.60 )- 
( [3.63 ) it follows that w kn _ k _ 2 are given by the matrix elements of the operators of dimension 
n + 3. This implies that "(terms involving u)f r _ l _ 2 \ r<n )" corresponds to matrix elements 
of operators involving total derivatives, which are given by linear combinations of terms 
~ (d.) n ~ r A kr _ k _ 2 (k = 0, 1, . . . , r — 2; r = 2, 3, . . . , n — 1). Thus we conclude, from ( |3.64|) 
and ( |3.66| ), that 

L-fc- 2 (°) = u(0)(-i < D .) n - k - 2 a u .gG u .(0)(u5 .) k d(0) + (total derivatives). (3.67) 

The operators Q kn _ k _ 2 (k — 0, 1, . . . , n — 2) have the same conformal spin and may mix 
with each other, although they do not mix with Of r _/_ 2 for n ^ r. To write down the 
corresponding RG equation, it is convenient to introduce combinations that are even and 
odd under the substitution k — ► n — k — 2: 

SL = 9 ^- fc - 2 1 Q "" fc ' 2 ' fc (k = 0,1,...,^), (3.68) 



where 



/ ■ . 



n 



- 1. (3.69) 



It is straightforward to see that S^. k and S~. k possess opposite "parity" under the G-parity 
transformation. Therefore, these two sets of the operators do not mix with each other. The 
RG equations are given by 

+ ^ (0; = g (^) w ^ (0; ^ (3 - 70) 



24 



where the one-loop anomalous dimension (r^~) is a [n/2]x[n/2] ([(n— l)/2]x[(n— 1)/2]) 
matrix, describing the mixing. Note that the number of independent operators and thus 
the size of the mixing matrix increase with conformal spin. 

In order to determine the anomalous dimension matrices in ( |3.70| ) we make contact 
with inclusive processes. Similarly to the discussion of the Wandzura-Wilczek part, we 
consider the case where the operators S^, k are flavour-diagonal and take the forward matrix 
element of (|3.70| ) in between nucleon states \N(P,S)). The total derivatives in (|3.67|) drop 
out, and ( |3.7U| ) reduces to the RG equations for 



(N{P, S)\4>(iD .) n a v -gG v .{iD ±(k^n-k- 2)\N(P, S)} } 



familiar from studies of the evolution of the n-th moment of the nucleon parton distributions 
e(x,fi 2 ) and h L (x,fi 2 ) defined in ( |2.14| ) and ( |2.13| ). It is straightforward to see that these 
operators for the upper and lower sign, which are even and odd under k n — k — 2, 
coincide exactly with the basis employed in renormalizing the parton distribution functions 
e(x,fi 2 ) |22[] and h L (x,fi 2 ) plf , respectively.^ Matching with these results, we obtain 



r 



-iT— 



k.l 



k.l 



-Yi 



k+2,l+2 



-x. 



k+2,l+2 



(k,P- 
(k,l 



0,1, 



0,1, 



i K r. 



(3.71) 
(3.72) 



where Yjj and Xij are the mixing matrices in the notation of [^] and pl| |, i.e. are given by 
Eqs. (3.12)-(3.16) of and Eqs. (3.14)-(3.16) of [0], respectively. By solving (|3~70l) , we 
obtain 



(3.73) 



and the matrix elements of S^. k are related to u>K n _&_ 2 ] an d uJk,n-k-2} °f (|3.49|) and ( |3.50| ) 



as (see (p4j) and flOSp ) 



Ip^fk^-k^} 



\n— k+1 



l^k\(n -k — 2)\ 



(0\SU0^ 2 )\p-(P,\)): 



(3.74) 



^3 Ct; {fc,n-fc-2} 



N2 



-1 



\n—k 



l^k\(n — k — 2)\ 



\S±(0;» 2 )\p-(P,\)). 



(3.75) 



Here the upper (lower) superscript should be understood for n — 2, 4, 6, . . . (n — 3, 5, 7, . . .) 
on the right-hand side. The results (g^ -(g?7|), combined with (|3732|) , (|3T33p , (|3T47| ), 
( |3.48|) and ( |3.59|) give the /i 2 -dependence of the three-particle contributions h^ 9 (u, /i 2 ) and 
h\" )9 (u,fx 2 ) (recall that (I P in ( P^7| ) and ( gig ) is given by 



)). The evolution of the 



5 For the case of the parton distribution /i 2 ), the relevant quark-antiquark-gluon operator has an 

additional 175 inbetween the quark fields. The evolution of the corresponding operator is not affected by 
the insertion of ^75 p5| . 
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three-particle twist 3 distribution amplitude T(a) of ( |3.35| ) is thus specified completely. For 
the first few moments we obtain from ( |3.71|) - (|3.75|) :fl 



if: 



fl(Q 2 ) = LW'^VOfZQfi) 



3>?, 0] ) (Q 2 ) = (/ 3 >g, 0] 



(3.76) 
(3.77) 



'3/ 3 >J 



3PH0.3] 
T T 
3p U [2,l] 



fT T 

. h P u [2,i] 



/41ri I 15/^ 29 /m Hr» \ 
V 15°^ ~ 12°^ 10°^ + T2°G J 



(3.78) 

corresponding to n = 3,4,5. Here Cq = N c , and cj^-n = 28/3 is substituted in ( ft.76| ) 
(see ( |3.51| )). Substitution of (|3.76| )- ([3.78| ) in ( |3.51| ) determines the scale dependence of H 9 
(n = 3,4,5). 

The results in ( |3.71p -( |3.78p illustrate a complicated mixing pattern characteristic for the 
higher twist operators. In particular, H 9 (Q 2 ) for n > 5 {h 9 n (Q 2 ) for n > 4) are not directly 
related to H^ji 2 ) (h 9 n (fi 2 )), in contrast to (|3.58|) for the twist 2 operator. 

There exist, however, two important limits, N c — > oo and n — > oo, where the three- 
particle coefficients H 9 and h 9 n obey a simple evolution equation. The mathematical reason 
for this simplification is the same as for the similar simplification observed for the nucleon 
parton distributions /i^(x,/i 2 ) and e(x,fi 2 ) in | 25| , [22] , [23]] . To show this, it is convenient to 
express H 9 and h 9 n directly by matrix elements of S^. k , by substituting ( |3.74[ ) and ( |3.75| ) 
into (|3~47| ) and (jEffi) . We find 



for n 



and 



2,4,6, 



n 




[n + l)\ 




N T 

n 


m p 


[n + l)\ 




N T 

n 


m p 


[n + l)\ 




N T 

n 


nip 


[n + l)\ 



E 

fc=0 



1 - 



2(k + 2)' 



EC 

fc=0 



n 



(0|^- fc (0;/. 2 )|p-(P,A)), (3.79) 



(0|5+ (0;/i 2 )|p-(P,A)), (3.80) 



£<0|S+ (0;/i 2 )|p-(P,A)), 



(3.81) 



Eh 

fe=o V 



2(k + 2)' 



n 



(0\S- k (0^ 2 )\p-(P,\)), (3.82) 



for n = 3, 5, 7, . . . Setting p 2 = Q 2 in fl3.79|) - ([3.82|) and substituting (|3.73|) into them, we 
would reproduce the complicated mixing discussed above. However, in the large N c limit, 



6 It is worth noting that the evolution of f3pwj k n ^ k ~2} coincides exactly with that of the coefficients, 
appearing in the conformal expansion of the twist 3 three-particle distribution amplitude of the pion (?[] . 
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that is, neglecting terms 0(1/N%) in the anomalous dimension matrices r„ of ( |3.T1 ) and 
( |3.72|) , the following exact relations have been derived p5| j2~2"f : 



k=0 



In 



in 



3,5,7,...), 



(3.83) 



and 



k 

where 



k=0 



El 1 -9^) ( r n + ) fe = (1-9^ 



< / 2(fc + 2) \ / j. 



7 T+ 



/ 2(l + 2)\ 
M I n + 2 



(n = 2,4,6,...), 



(71 = 3,4,5,6,...), 



(3.84) 



(3.85) 



^ = 2NA^|J(n + l)+ lE -- A 



A 2(n + l) 



(3.86) 



2N C U(n + 1) + 7s - 7 + 



2(n+l) 



As a consequence of these relations, we obtain 

(Q 2 ) = rf ± ' b (f?h° n ) (^ 2 ), 



(3.87) 



(3.88) 
(3.89) 



where the upper (lower) superscript should be understood for n — 2, 4, 6, . . . (n = 3, 5, 7, . . .) 
on the right-hand side. Therefore, in the large N c limit, and /i^ obey simple DGLAP-type 
evolution equations similarly to the twist 2 case Q3.58 ); they are governed by the anomalous 
dimensions given in analytic form in (3.86) and ( 3.87| ). We note that these anomalous 
dimensions correspond to the lowest eigenvalues of the mixing matrices [|~ 
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The phenomenon leading to ( |3.88| ) and (|3.89|) can be stated as decoupling of the three- 
particle operators, which correspond to the higher eigenvalues of the mixing matrix, from 
the RG equation. The same decoupling is observed at large n for arbitrary values of iV c 
23fl . In this case, we obtain ( |3 .83| ) — (|3. 85| ) with the anomalous dimensions ( |3.86| ) and 



( ET87D shifted by 



7 T± 

In 



7^ + (4C F - 2N C 



In n + ^ E - - 



(3.90) 



With this modification of the anomalous dimensions, the results (|3.88|) and (|3.89|) are valid 
to O ((1/Nc) \n(n)/n) accuracy. 

These simplifications provide useful insight both into the model-building of the distri- 
bution amplitudes, and a convenient description of their scale-dependence: In the large 
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iV c -limit, each conformal partial wave of the three-particle contributions is described by a 
single nonperturbative parameter, as demonstrated in ( |3.88|) and ( |3.89| ). This is remarkable 
because in general each conformal spin involves many independent nonperturbative matrix 
elements (see ( |3.79| )-( p.82j )). This point can be made stronger with the full account for 
effects subleading in N c but for large conformal spin j, as shown in ( |3.90 ). Furthermore, 
fl3.90| ) proves the conjecture made in [ID, |7| that the lowest anomalous dimension of the 
twist 3 three-particle operators is increasing as ~ lnj, similarly to the twist 2 case (|3.57| ). 
This ensures convergence of the series in the Appell polynomials at least for high energy 
scales. Combined with the fact that the distribution amplitudes can be resolved order by 
order in the conformal spin, the truncation of the conformal expansion at some low order 
provides a useful and consistent approximation of the full amplitude. 

To summarize, we have worked out the scale dependence of chiral-odd twist 3 distri- 
bution amplitudes h^\u, jj 2 ) and h\f (u, /i 2 ) in the leading logarithmic approximation. In 

the two limits, N c — > oo and n — > oo, the evolution of h^\u, /i 2 ) and h^\u, /i 2 ) is dras- 
tically simplified and reduces to a DGLAP-type equation. The discussion in this section 
completes the results for the chiral-odd distribution amplitudes up to twist 3, which can be 
predicted based on the QCD constraints from equations of motion, conformal invariance, 
and renormalization group invariance. 



4 Chiral-even Distribution Amplitudes 

The analysis in the previous section can directly be extended to the chiral-even distribution 
amplitudes. In Sec. 4.1 we derive the constraint equations for g±(u) and g±{u) imposed 
by the QCD equations of motion, and identify the contribution to g±\u) and g±\u) from 
the twist 2 distribution amplitudes (p\\(u) and 4>±(u) and the three-particle distribution 
amplitudes V(a) and A(a). In Sec. 4.2, we study the conformal expansion for 4>\\(u), g± (u) 
and g±\u). In Sec. 4.3, we work out the renormalization of g±\u) and g±\u), utilizing the 
fact that the conformal symmetry is preserved at one-loop level. Our presentation in this 
section will be brief, since the methods and the results are in parallel with those in Sec. 3. 



4.1 Equations of motion 

To derive the constraint relations among the chiral-even distribution amplitudes we again use 
operator identities (to twist 3 accuracy) for the nonlocal operators in (|2.19|) and (|2.20| ) JTOf: 

r 1 d _ 

ufejjulx, — x)d(— x) — / dt — — u(tx)p\tx, — tx]d(— tx) 

J CyXjj^ 

— / dtt / dvu(tx)[tx, vx\gG lJbU {vx)x l '^ r y^[vx, — tx\d{— tx) 
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fl rt 

— i I dt dvvu(tx)[tx,vx]gG^ l/ (vx)x I/ [vx,—tx]^d(—tx) 



— it/ivap / dttx v d a u(tx)'y l3 'y5[tx J —tx]d(—tx) 



and 

m(^)7m75[^, -x]d(-x) 



+ (m u — m d )x u I dttu(tx)a l/tM [tx, —tx]d(—tx), 



1 d 

dt — — u{tx)^^\tx, —tx]d(—tx) 
ox^ 



(4.1) 



dtt J dv u(tx)[tx, vx]gG flu (vx)x u qt[vx, — tx\d{— tx) 

— if dt f dvvu(tx)[tx,vx]gGnv(vx)x u [vx,—tx]qt'y5d(—tx) 

Jo J-t 

— itfiva/3 / dttx u d a u(tx)^[tx, — tx]d(— tx) 



+ {m u + m d )x v / dttu(tx)a I/ll 'j 5 [tx } —tx]d(—tx), 



(4.2) 

where d a is the total derivative defined in (|3.3j), and the terms proportional to quark masses 
originate from the use of QCD equation of motion. By sandwiching these equations between 
the vacuum and the p meson state, and taking the light-cone limit x — > z, one obtains the 
following relations among the distribution amplitudes: 



due^g^iu) = I' dt j~ du^ pz (j)\\{u) - (£(pz) 2 I t 2 dt /" dvA(v,tpz) 



i ,i 

dt / 
ii Jo 

l 



- Cl(P z ) 2 L dtt2 I dvvV(v,tpz) - -(pz) 2 (1 - 5+) I dtt 2 I du^ pz gf{u) 



id-ipz) f 1 dtt f 1 due^ pz ^ ± (u), 
Jo Jo 



1 

-1 
1 



D it£pz (a). 



(4.3) 



and 



-(l-5+) f 1 du e&'gf (u) = ['dtt t du e^gV (u) 
2 v ' Jo ~~ Jo Jo 



+ iQ(pz) dtt 2 dvV(v,tpz) +i(t{pz) dtt 2 dvvA(v,tpz) 



1 rl 



dtt / du(? tipz (j) ± {u), 
o Jo 



(4.4) 
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where V(v,tpz) and A(v,tpz) are defined in Q2.27] ) and we introduced the notations 



fT 2 fT I rV,A 

r _ J P x Jp m u =t rn d V A j 3p 

°± = - j t 2 -6 ± = -f- ' C 3 ; = ~r—- ( 4 -5) 

Jp Jp rn p fpf^p 



In order to solve these equations, we expand them in powers of pz and transform them into 
relations among the moments of distribution amplitudes: 



n(n — r 



,71 



+ 1)MW = Mil + -L— -1 (l - 6+) M { n % + C 3 >(n - 1) dvAn- 2 (v 



-l 



+ C 3 >(« - 1) JdvvV n - 2 (v) - LnM^, (4.6) 
-l 

and 

1 i i 

-(n + 2) (l - 5+) = + (V n JdvV n -i(v) + C 3 > jdvvA n ^(v) - 5+M^ 



-l 



(4.7) 

where V n (t>) and A n {v) are defined similarly to (|3.10| ) from V(a) and ^4(a), and we intro- 
duced the shorthand notations M^'^ = ^du^ n . From these equations, one obtains 
a recurrence relation for M!£> as 

*) 



(l - 6+) ((n + 2)(n + l)Mjfi - n(n - 1)m£> 

= 2Mj + 2C 3 y p [n(n + l)V n ^(v) + n(n - l)vV n - 2 (v)] 



+ 2C 3 A P / dv [n{n + l)vA n ^{v) + n(n - i)An- 2 (v)] 



53/3 



-1 



2(n + l^+M 1 / - 2n5_M I /_ 1 . (4.8) 



This equation is similar to ( |3.14 ) for (u) , and can be cast into the form of a differential 
equation as 

(1-6+) uu(g^)"(u) = -*(«), (4.9) 

where 

*(«) = 20|| (u) + 5+^'Au) + 6^' ± (u) 
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d U U I / 

+ 2C,L -T- da d da u - 

F du J J 1 — ad — a u \ 



o o 



d u u I / 

+ 2(£ ) — da d da u - \a d 

F du J J 1 — ad — a u \ 



oca 3 h «« -j — 

dad da u 



CM, ,, 







dad u da 



From this equation, one immediately obtains the solution for g± (u) as 



v(s) 



,A(a). (4.10) 



7 1 } 1 

(l - <5+) #[ a) (w) = u Jdv- V(v) +u Jdv- $(v). 



(4.11) 



Combining this result with ( |4.7| ), one obtains the solution for g±\u) as 

u 1 



^(«) 



(if - V(v) + dv - V(v 
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+ C A /da /da 1 f ^ + ^ 

3p i i u 1 — ad — a M \dad da Ui 

d }. } . V(a) 



A(a) 







(4.12) 



Eqs. ( j4.ll ) and ( |4.12 ) again allow the decomposition of g±'(u) and g±\u) into several terms 
according to the various source terms: 



#W = 



(4.13) 
(4.14) 



where gfj" (u) and g^ WW (u) denotes the contribution from the twist 2 distribution am- 
plitudes (Wandzura-Wilczek part), g± (u) and g^ g (u) are the contribution from the three- 
particle distribution amplitudes V and A. In particular, we get 



(v)WW , x 



g l f >ww (u) = 2u I dv -4>\\(v) + 2u I dv —<fi\\(v). 



v 



1 , 1 



dv -<f>\\(v) + / dv —<j>\\[v) 



v 



1 , 1 



(4.15) 
(4.16) 
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4.2 Conformal expansion 

The conformal expansion for the chiral-even distribution amplitude can be performed sim- 
ilarly to the case for the chiral-odd ones in Sec. 3.2. In the following, we restrict ourselves 
to the massless case. 

For completeness, we start with the expansion for the twist 2 distribution amplitude 
4>\\{u). It reads 

oo 

0n(«) = 6ufl£4c2 /2 (O. (4.17) 

n=0 

where C^' 2 (£) is the Gegenbauer polynomial, and each term corresponds to the conformal 
spin n + 2. We note a| = 1 due to the normalization condition ( |2.23fl . Because of the 
G-parity invariance of the p meson distribution amplitude (likewise for other mesons such 
as to, 4>), it follows that a| = for n — 1, 3, 5, . . . in ( 4.17|) . In Sec. 5, however, we shall 
treat an application to the K* distribution amplitude with explicit SU(3) breaking due to 
the quark masses. With this in mind, we keep all terms a| in (|4.17|) . For the same reason, 
we shall work out the conformal expansion for g*f\ V(a) and .4(a) by keeping both 
G-parity invariant and G-parity violating terms in the following. 

To carry out the conformal expansion for g±^(u) and g± (u), we again introduce a set 
of auxiliary amplitudes g^{u) and g^(u) [^(J defined by 

(Q\u(z) r ^[z, -z}d{-z)\p-{P, A)) = -f p m p e^ f du^ z g^(u) } (4.18) 

J 

(0\u(z)^j.[z, -z]d(-z)\p-(P, A)) = -f p m p e% f du^g^(u), (4.19) 

J 

which are related to g± (u) and g^ (u) as 

gHu)=gi\ U ) + ~gi\u) t (4.20) 

g^{u) = gt\u)-- A ^-gt\u). (4.21) 
4 du 

The conformal expansion for g^(u) and g^'(u) is given by 

oo 

» T1 («)=2HE»W )) (0, (4-22) 

ra=0 
oo 

^ t («)=2«i;^ ' 1) (0- (4-23) 

n=0 
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Substituting these expansions in ( |4.2(J| ) and ( ^.21| ), one obtains 



9±W 



E (Gn - g„_i) cy 2 (o + e (^-^-i)cy 2 (o, 

n=0,2,4,... n=l,3,5,... 



(4.24) 



( a ) / ^ 
9±(u) 



E , G r i ^To^ n /2 (fl+ E TT^T^r^ gyCO 1(4-25) 



v n=0,2,4,.. 



(n + l)(n + 2) 



n=l,3,5,. 



(n + l)(n + 2) 



where G n and g n represent, respectively, G-parity invariant and G-parity violating compo- 
nents in the expansion defined by 



G„ = 



g n + (-l) n gi 



IT 



- (-l)^IT 



(4.26) 



for n = 0, 1, 2, . . ., and G_i = g_i = is implied. Note the difference between {G n , g n } and 
{H n , h n } (see ( |3.34|) ) owing to the chiral-even or -odd nature of the distribution amplitudes. 

The conformal expansion for the three-particle distribution amplitudes V(a) and A(a) 
can be written down similarly to (|3.35| ): 



V{a d , a u , l-a d -a u ) = 360a d a u (l -a d - a u ) 2 E ^k,i J k,i{ a d, at u ), 



k,l=0 



A(a d , a u , 1 - a d - a u ) = 360a d a u (l - a d - a u ) 2 E ^k,iJk,i( a d, (4.27) 

k,l=0 



,v 



V 



The G-parity invariance of the three-particle distribution amplitudes leads to LU kl = —uj l<k 
and ui kl = ujf k . As was stated at the begining of this subsection, we shall not assume 
this symmetry in the following. We also note the normalization condition in (|2.29|) gives 



UJ 



[0,1] 



28/3 and Uq = 1. The conformal spin for each term in the above expansion is 
equal to j = + / + 7/2, and the preservation of the conformal invariance at one- loop level 
prevents mixing among the contribution with different n = k + I. 

Our next task is to identify the twist 2 (Wandzura-Wilczek) and the three-particle con- 
tibutions to G n and g n . We decompose 



G n — G 



WW 



+ GI- 



9n = g 



WW 



(4.28) 



and consider the Wandzura-Wilczek contribution first. Substituting (|4.15| ) and Q4.16D into 
flCTD and fl£2Ip and using the formulae ( [7021) , fl/Hp, (|A~^ ) and flO|) , we obtain the 
Wandzura-Wilczek contribution for g^(u) and g^(u) as 



WW, 



u 



i dv m 
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oo 

n=0 



3(n + 2) 
2(2n + 3) 



3(^+1) p(i,o) r ^ 
2(2n + 3) n+lW 



(4.29) 



at 



WW, 





oo 

2w Yl a l 

n=0 



2(2n + 3) 

These give rise to G^ w and g^ w in ( |4.2g| ) as 



2(2n + 3) 



(4.30) 



r ww = 3(n + 2) || 
n ~2(2n + 3) n ' 



-3n 



2(2n + 1) 



4 n— 1) 



G 



n 



-3n 



2(2n + 1) 



3(n + 2) || 
2(2n + 3) an ' 



(n = 0,2,4 



(n = l,3,5...) 



(4.31) 



For even n, we find that a\ which corresponds to the conformal spin n + 2 in the expansion 
for the twist 2 distribution amplitude gives rise to G^ w and G™7{ which corresponds to the 
conformal spin n + 3/2 and n + 5/2, respectively. Likewise for odd n, a| gives rise to g^ w 
and • This is the same pattern as observed for the chiral-odd distribution amplitudes 
and is discussed in App. [B|. 

From the solution for g± (u) and g±\u) in (|4.11|) and (|4.12 ) we can identify the three- 
particle contribution to the auxiliary amplitude as 



g U9 (u) = Cl { l l dv -H(v) + M(u)\ + (i p { f dv -L(v) + 

wit V J {Ju V 

d [£ dv ~H{v) + M(«)} + Ci {[ dv ±L{v) + N{ 



g U9 (u 



where the functions H(u), L(u), M(u), N(u) are defined as 

d d 



(4.32) 
(4.33) 



H u] 



d r u 

du 
d 



ru I 

dad / da u — 
Jo a„ 



a d - h a u 

dad 



du 



d 



d 



L{u) = — / da d / da u — a d - a u 



a B 



da r , 



d 

'da,. 



Via), 
A(a) 



du 



M(u) = — / dad / da u — V(a 



a n 



ru ru | 

N(u) = / da d / da u — 
J o J o a a 



d d 
+ 



da d da u 



A(a) 



(4.34) 
(4.35) 
(4.36) 
(4.37) 
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with a g = 1 — ad — ot u . The calculation of ( 4.32| ) and ( |4.33| ) can be carried out similarly 
to the chiral-odd distribution amplitudes in Sec. 3, using the explicit form for V(a) and 
A(a) in (|427D . First H(u), L(u), M(u) and N(u) in (|Q1 - (^37D can be easily obtained by 
using ( |A.15| ) — ( |A.18| ). The integration of H(u) and L(u) in ( f4 . 3 2| ) and ( |4.33| ) can be done by 

to the results of those integrations and N(u), and the 



( [TOD . Application of ( \AJ§ and 
use of ( |A.5[ ) and ( |A.6[ ) for M(u) yields the anticipated form for the result. One eventually 
obtains 



x{n-(2k + 2-n)} (-C 3 >^_ fe _ 2 + C 3 >t 

oo n— 2 

90(i + OE^ 0,1) (OE 



(4.38) 



n ~ 2 £;!(n- A; - 2)!(-l) fc 



n=2 



k=0 



(n + l)!(n + l) 



x {n + ( 2 fc + 2 - n)} + C^_ fc _ 2 ) 



(4.39) 

where we kept the form (2k + 2 — n) in {. . .} which is anti-symmetric under the interchange 
k n — k — 2. From this equation we get G 9 = G{ = g% = gf = 0, and 



Gt 



^ k\(n - k - 2)\(-l) n+k 

90 — TT i ( - A' + - - /' it, ; 



A _ (1 (n+l)!(n + l) 

(n = 2,4,6...) 
«Z?£;!( n _£;_2)!(-l)"+* 

90 7 ; rT7 ; 

to (n+l)!(n + l) 
(n = 3,5,7...) 

^ 2 A;!(n-A;-2)!(-lf+ fc + 1 

90 > ; --— ; 

fro (n + l)!(n + l) 

(n = 2,4,6,...) 
n ~ 2 £;!(n-A;-2)!(-l) n + fc + 1 



{(2fc + 2-n)C ; 



3p u '[fc,n-A;-2] 



21 + <3p^ 



3p uy {fc,n-fc-2}/5 



3p U; [fc,n-A:-2] 



+ (2A; + 2-n)C 3 >f fc , n _ fc _ 2 }} J 



{<• 



3p U Jk,n-k-2} + (2A: + 2 — ^)Cst Ct; [fc,n-A:-2] 



(4.40) 



}■ 



90 £ 

fe=0 



(n + l)!(n + l) 
(n= 3,5,7,...) 



{(2fc + 2 - n)C]^{ifc in _ fc _2} + < 3 Vffc,n-fe-2]}' 



(4.41) 



where we introduced the anti- symmetric and symmetric components of uj^'f defined by 



v,a 



V,A V,A 



k.l 



I,k 



(4.42) 
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One finds here again that G n and g n which correspond to the conformal spin j = 3/2 + n 



receives contributions from the coefficients uj^'f with a fixed k + I = n — 2 which have the 



same conformal spin j = 3/2 + n. The first terms in fl4.40|) and ( [4.41D read 



G% = IOC3X0 = 10C 3 A P , 



Gl - - (— 2C^[o,2] + Kipufop} ~ Ct P u {i,i} ) ■ 



gl = jj (2C 3 >£ 2] - 4C 3 >f , 2} + C3>[i,i } ) , • • • , (4-44) 

where we have used the normalization condition for A(a) and V(a), cu^ = 1 and w?oi] = 

28/3. We also note that ( |4.31| ) for n = 0, 1 together with the conditions ajj = 1, Gq = G\ = 
determines the first two coefficients of G n as Go — 1? G± = —1/2, which gives consistent 
normalization ^dug^\u) = 1 and J^dug^f'^) = 1 through fl4.24|) and ( |4.25|) . We finally 



remind that the G-parity invariance of the three-particle distribution amplitudes imposes 
V{k,i} = and rf kA = 0, leading to g% = and G 9 n = gH°. 

4.3 Renormalization and scale-dependence 

In this subsection, we discuss the renormalization of the chiral-even distribution ampli- 
tudes, 4>\\(u, /i 2 ), g±(u, /i 2 ) and g±\u, /i 2 ), utilizing the conformal expansions derived in the 
previous subsection. 

For completeness we start our discussion with the renormalization of the twist 2 dis- 
tribution amplitude <f)\\(u). From ( ^4.1| ) and the orthogonality relations of the Gegenbauer 
polynomials [[HJ, one obtains 



= T( 2 f " | 3 | ! l duCT{0Un,H 2 ). (4.45) 
3(n + l)(n + 2) jo 



Using ( |2.19| ) in the right-hand side of (|4.45 ), we can express a\ in terms of the conformal 
operator: 

a l^ 2 ) = T} 7 ZT^rZWlX) V vr<0|^ll(0;^ 2 )|p-(^, A)> (4.46) 

6jpTn p {n + lj(n + 2 l )(e (A < , • z)(j> • z) n 
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with 



Ql(x;^) = {id.Yuix^.CT ( |l ] d(x), (4.47) 



where the local operator in the right-hand side is renormalized at /z 2 , D= D — D , and <9 M is 
the total derivative ( |3.3| ). fi|[(x; li 2 ) is the conformal operator with conformal spin j = n + 2 
TEfl . The scale dependence of a\{jj?) is well known |], 0: 



ai(Q 2 ) = L^ b ai(/i 2 ), (4.48) 

where L = a s {Q 2 ) / a s {ii 2 ), b = (lliV c — 2iV/)/3 and the anomalous dimension 7JI for the 
conformal operator fi^ is given by 

7 « = 4C F („(„ + 2)+ 7E -|- 2(n+i ; (n + 2) ) • (4.49) 

For n = 0, f)JJ * s reduced to a conserved vector current, and hence its anomalous dimension 
vanishes to all orders. Combined with the normalization condition for (j)\\(u), ajj = 1, this is 
consistent with the fact that f p is scale independent. We thus omitted f p in both sides of 
flPgD (compare with (^58|) and (|3T59|) ). 



Next we proceed to discuss the renormalization of the twist 3 distribution amplitudes 
g±\u, jj 2 ) and g± (u, /j 2 ). As we saw in the previous subsections, they receive contributions 
from the twist 2 distribution amplitude (Wandzura-Wilczek parts), the three-particle distri- 
bution amplitudes and the terms proportional to the quark masses (see ( [113| ) — ( (4.16| )). The 
scale-dependence of 0||(w,/i 2 ) discussed above completely determines that of g^ ww (u, /i 2 ) 
and g ( f )WW '(u, /1 2 ) through the relations ( ggg ), ( ggg) and ( ggT]) . 



To understand the scale dependence of the three-particle distribution amplitudes V(a) 
and A(a), one needs to express uj^'f in ( [4.27] ) in terms of the local conformal operators. 



Owing to the orthogonality relation (|A.13|) of the Appell polynomials with different k + I, 



lo^\ can be expressed in terms of the matrix elements of the conformal operators with a 
definite conformal spin j = k + I + 7/2. As we saw in Sec. 3, it suffices to know the form of 
the conformal operators up to total derivatives for the renormalization. For this purpose, 
we recall from ( |2.24| ) and ( |2.25| ) 



(0\u(tz)j.[tz, vz]gG±.(vz)[vz, wz]d(wz)\p (P, A)) 

= ie^p^lip ■ z)fl J Vae- l ^ taM ^V{a), (4.50) 
(0\u(tz) r y.i r y 5 [tz, vz\gG±.(vz)[vz, wz\d(wz)\p~ (P, A)) 

= ie^vP^lip ■ z)fi p J Vae- ip < taMd+va ^A(a). (4.51) 
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To obtain the actual form of the conformal operators, we apply d n ~ 2 / dt n ~ k ~ 2 dw k on both 
sides of (|4.50|) and ( |4.51|) and set t = w = v = 0. Using the integral formula for the Appell 
polynomial (|A.14j) , one obtains 



(0\u(0)(iD . 



,n-k-2 



7 .gG ± .(0)(iD.) k d(0)\p(P,\)) 



\n-l 



V 



360(-l) n A;!(n - k - 2)! 
2"+ 1 (n + l)(2n + 1)!! 



+ (terms with ujY r _ l _ 2 \r<n) 



(4.52) 



(0|ti(0)(< D T^l-^gGUm D .) d(0)\p(P, A)) 



A 360(-l) n k\(n-k-2)\ 
M-fc-2 2"+ 1 (n + l)(2n + l)!! 



+ (terms with co , / 4 r _;_ 2 |r<n) 



(4.53) 



It is easy to see by induction that wYr^-i-2\r<n m these equations are the matrix elements 
of the total derivatives of the lower conformal operators. Therefore we can identify the 
corresponding conformal operators for (jJ['n-k-2 as 

(-!)*#, 



(/ 3 Xn-fc- 2 ) (^) " gm{n _ k _ 2)[ 



(0|e^ n _ fe _ 2 (0;^)|p-(P,A)), 

(o|e^_ fc _ 2 (o ;yU 2 )| P -(p,A)), 



(4.54) 
(4.55) 



90A;!(n - Jfe - 2)! 

where the conformal operators are now obtained up to total derivatives as 

9^ n _ fc _ 2 (0) = u(0)(-i D .) n - fe - 2 7.#<5±.(0)(i D .) k d(0) + (total derivatives), (4.56) 

@fcn-fc- 2 (0) = u(0)(-i 5 .) n ~*~ 2 7.</Ga-(0)(i B 0^75^(0) + (total derivatives), (4.57) 



and we introduced for convenience a dimensionless and scale independent normalization 

2 n -\2n + l)!!(n + 1) 



constant iV n as 



(4.58) 



From Eqs. (|4.54j) and Q4.55|) , we can obtain the scale-dependence of (/^<^ n _fc_ 2 )(/i 2 ) and 
b y working out the renormalization of {6^ n _ fe _ 2 , ©fe jn _ fc _ 2 } (k = 0, . . . , n- 
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v± 



@fc^n-fc-2 i ©n-fc-2,fc 



0,1,..., with defined in ( |3.69p ) 
0, 1, ■ ■ ■ , {i^J , R££} and {i^, i^Tj, respec- 
l) n+1 and (— l) n and thus they do not mix with each other under 



2).[] If we define R v n 

and R n j. = Qk,n-k-2 T ©n-fc-2,Ai 

tively, have G-parity 
renormalization. 

By inserting (|4.54 ) and ( |4.55| ) into ( |4.40| ) and ( |4.41| ) and recalling the definition of C^ p ,A 
from (|4.5| ), one can express the contribution from the three-particle distribution amplitude 
to the two-particle distribution amplitudes g^'^ in terms of the conformal operators: 



N n 


f p m p (n + 1)! 


(n + 1) 


-N n 




f p m p (n + 1)! 


(n + 1) 


-N n 




f p m p {n + 1)! 


(n + 1) 






f p m p (n + 1)! 


(n + 1) 



n-2 



-(0|^(fc + l)it:- fe (0;// 2 )|p-(P,A)} ) 



fc=0 
n-2 



(0|£(n-*-l)J^(0;,i a )b-(P,A)), 

k=0 



n 



n 



2,4,6, 



3,5,7,...) 



(4.59) 



n-2 



(0\J2(n-k-l)Rl k (0;fi 2 )\p-(P,X)), 

k=0 
n-2 

(0\Y,(k + i)R- k (0;» 2 )\ P -(P,\)), 



in 



in 



2,4,6, 



3,5,7,...) 



fc=0 



(4.60) 



where the operators R^ k (k 



0, . . . , n — 2) are defined as 



-2 ^ ®n-/c-2,fc + ®fc,n-fc-2 + ® 



n—k—2,k' 



(k = 0,1,..., n-2) (4.61) 



Here we note that {-R^} and {i? n fc } have opposite parity under G-parity transformations 
and constitute another operator basis which is equivalent to {R^Xt R-nX) an< ^ i-^n k > ^nfc}) 
respectively, either of which has the same number of independent operators + k~ + 2 = 
n — 1. If we define the anomalous dimension matrix for {R^k} as T^, the scale dependence 
of R^ k is given by 



n-2 



Rik(o;Q 2 ) = Y.{L r " /b ) kl RnM» 2 )- 



(4.62) 



1=0 



The renormalization for {-R^} can be conveniently worked out by considering the for- 
ward matrix elements with respect to a spin 1/2 target (say a quark) as was done in 
Sec. 3. In this case, the contribution from the total derivative terms disappear and its 
matrix element (R^ k ) is reduced to (ip(iD.) n ~ k ~ 2 'y.gG±.(iD.) k ifj ± (k — > n — k — 2)) + 
{ip(iD .) fc 7 .gG ±.(iD .) n ~ k ~ 2 i^il) + (A; — > n — k — 2)). By now, renormalization of {R nk } 
has been solved by several different approaches in the context of the Q 2 evolution of the 
transverse spin structure function g2(x, Q 2 ) [110, ^j], In particular, {R^k} is exactly the 



7 since er 



and 0jf„_fe_ 2 have the same conformal spin, they generally mix with each other under 



renormalization, even though they originate from different distribution amplitudes 
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operator basis used in [E7 for the renormalization of g2{x,Q 2 ). From the result in P7| , r+ 
in ( |4.62| ) is obtained as 



(rj) 



(fc,/ = 0,l,...,n-2) 



(4.63) 



where X^ 1 is given in Eqs. (15)— (17) of J27| with n — > n + 1. The explicit form of T~ 
is not available in the literature, but can be obtained from the kernel given in [TIJ. (Since 
R~ k does not contribute to the deep- inelastic scattering which is the charge conjugation 
even process, it has not been receiving attention up to now.) 

By using the basis {R„k} an d {R-n fc}> (|4-54|) and ( [4.55| ) can be rewritten as 



(fl 



3p a; [fc,n-fc-2] ="= ^3p u '{fc,n-A:-2} 



:-D k N„ 



180fc!(ra-fc-2)! 



(0\Rl k (0;^)\p-(P,X)), (4.64) 



for the G-parity invariant components and 



A, .A 



3p U \k,n-k-2} ~F f£p U \k,n-k-2] 



'-l) k N„ 



180fc!(ra-fc-2)! 



(0K fe (0;/i 2 )|p-(P,A)), (4.65) 



for the G-parity violating components. In ( ^4.64| ) and ( [4.65| ), upper and lower signs corre- 
spond to n = 2, 4, 6, . . . and n = 3, 5, 7, . . ., respectively. For illustration, we give here the 
explicit form of the scale dependence of the left hand side of ft4.64|) for n = 2,3,4 from 
( jOg) : 



fi P {Q 2 ) = L T ^ b f 3 A p (v 2 ), r 2 - = --w 



-CV + 3C, 



Gi 



(4.66) 



3 i3p i; 

28 fV , f 

3 i3p ' J; 



A .A 
3p W {0,l} 



A, ,A 
3p UJ {0,l} 



3 i3p J; 

28fV i f A 

3 ^3p i: 



^{0,1} 



3p^{0,l} 



IC7. 



G 



\Cg rCf + 4Cg 



(4.67) 



a uj a 



3p u/ {0,2} 



/ 2/ 3 >[o,2] + 2/. 

-/3p W {l,l} 

V -2/ 3 > [0 ,2] + 2/ 3 > 



L 1 * 



3p u/ {0,2} . 



fA, ,A 

~J3p UJ n,i} 



3p^{l 

V — 2/ 3p a;[o ) 2] + 2f 3p u {02} 



G 



31 j_ 25^ 

-c F + la 



12 



a 



23 , 



2"Oi? — ^Or 




(4.66 



8 The result in Eqs. (15)-(17) of is correct only for odd n in their convention. For general n, the term 
(-l)V(ri - proportional to 2C F - C G in Eq. (16) should read - /). See @ for the detail. 
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where we have used Uq = 1 and ojY q ^ = 28/3, and the anomalous dimensions T 2 , T 4 can 



be obtained from the kernel given in [10|. [] For the complete discussion on T n , see ||29|| . 

As one can see from this illustration, the u 2 evolution of {uZ n _ k _ 2 y u {k,n-k-2}} anc ^ 
{^{A; n _fc_2}! u \k n-k-2]\ becomes extremely complicated for general n. However, as was the 
case for the chiral-odd distribution amplitude in Sec. 3, the /^-dependence of the three- 
particle contribution to the twist 3 two-particle distribution amplitudes, i.e. G 9 n (n 2 ) and 
g^[(u 2 ), is greatly simplified in the large N c limit. It has been shown in [^] that the 



combination of R^ k in ( |4.59| ) and ( |4.6(J| ) renormalize multiplicatively at N c — > oo to 0(1 /N ( 



'2\ 



accuracy. We thus get in this limit 

Gi(Q 2 ) = L^ b G° n (f, 2 ), 

gi{Q 2 ) = L^ b g^ 2 ), (4.69) 
with a common anomalous dimension 

In = 2N C + 1) + lE - I + ^Vl) ) • ( 4 ' 7 °) 

This 7 n has been shown to be the lowest eigenvalue of ||26||. This reduction to the simple 



evolution equation ( |4.69| ) is equivalent to the fact that the coefficients of i?„ k in (|4.59|) and 



Q4.60| ) constitute the left eigenvector of r„ with the eigenvalue 7„ in this limit: 



n-2 



+ (r-) ki = (i + i) lr , 



k=0 

-2 



J2(n-k-l) (r+) fci = (n - / - 1) 7 „, (4.71) 

k=0 

which can be compared with (|3.83|) — ( |3.85| ) for the chiral-odd case. The renormalization 
of the flavour singlet part of g±'(u, u 2 ) (for u> and <p mesons) is complicated by addi- 
tional mixing with the purely gluonic twist 3 distribution amplitudes. For this mixing, no 
simplification occurs in the N c — > oo limit. 

As was discussed for the chiral-odd distribution amplitudes in Sec. 3, another simplifica- 
tion of the evolution equation for g^' a \u, fi 2 ) occures at large n with arbitrary N c . In this 
limit, the scale dependence of G^(/i 2 ) and g^(/i 2 ) is described by the same equation ( |4.69|) 
with a slightly shifted anomalous dimension 

In -> In + (4CV - 2N C ) (inn + lE - £) . (4.72) 



Combined with this simplification at n —>■ oo, the result at large N c in Q4.6S ) is valid to 
0((1/N 2 )\n {n)/n) accuracy as was the chiral-odd case in Sec. 3. 

9 Eq. (6.2) in for the kernel contains a misprint: the delta function in the last line should be replaced 
by S(a — u). 



41 



To summarize this section, we have solved the renormalization of the nonsinglet chiral- 
even twist 3 distribution amplitudes, g±\u) and g± (u). We found that in the limits N c — > 
oo and n — > oo, the scale-dependence of the three-particle contribution to g± (u) and g±\u) 
is described by a simple DGLAP type evolution equation similar to that for the twist 2 
distribution amplitude. Combined with the results of the previous section, this means that 
this simplification for the scale-dependence is universal for all twist 3 nonsinglet distribution 
amplitudes. 



5 Models for Distribution Amplitudes 

In this section we present explicit models for the two-particle distribution amplitudes of 
twist 3, taking into account contributions of the first few conformal partial waves. The 
main observation and important point is that the QCD equations of motion are satisfied 
order by order in the conformal expansion, which guarantees the consistency of the ap- 
proximation. Our approximation thus introduces a minimum number of nonperturbative 
parameters describing matrix elements of certain local operators between the vacuum and 



the meson state, which we estimate using QCD sum rules |5D, [3]]. More sophisticated mod- 
els can be constructed in a systematic way by adding contributions of higher conformal 
partial waves when estimates of the relevant nonperturbative matrix elements will become 
available. 

Our approach involves the implicit assumption that the conformal partial wave expansion 
is well-convergent. This can be justified rigorously at large scales, since the anomalous 
dimensions of twist 2 and twist 3 operators increase logarithmically with the conformal spin 
j, but is nontrivial at relatively low scales of order /x ~ (1 — 2) GeV which we choose as 
reference scale. We believe that this assumption is natural and in fact necessary for any 
model of distribution amplitudes at scales where they evolve perturbatively; the last word, 
however, has to come from experiment. 

Since orthogonal polynomials of high orders are rapidly oscillating functions, a trun- 
cated expansion in conformal partial waves almost necessarily is oscillatory as well. Such 
a behaviour is clearly unphysical, but does not constitute a real problem because physical 
observables are given by convolution integrals of distribution amplitudes with smooth coef- 
ficient functions. A classical example for this feature is the 77*-meson form factor which is 
governed by the quantity 

/ du-(p(u) ~ Vaj, 
J u 

where the coefficients a, are exactly the "reduced matrix elements" in the conformal expan- 
sion. The oscillating terms are averaged over and strongly suppressed. Stated otherwise: 
models of distribution amplitudes should generally be understood as distributions (in the 
mathematical sense). 
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V 


p ± K*± <f) 


/y[MeV] 


198 ±7 226 ±28 254 ± 3 



Table 3: Experimental values of couplings to the vector current [|31j] . 



V 


p ± K*± <p 


/£(lGeV)[MeV] 


160 ±10 185 ±10 215 ±15 


aj(l GeV) 
4(1 GeV) 
4(1 GeV) 
4(1 GeV) 


0.19 ±0.05 
0.18 ±0.10 0.06 ±0.06 0±0.1 
0.20 ±0.05 
0.2 ±0.1 0.04 ±0.04 0±0.1 



Table 4: The tensor couplings and lowest Gegenbauer moments of vector mesons from QCD 
sum rules, see App. 0. 



5.1 Leading twist distributions 



The twist 2 distribution amplitudes of vector mesons have received much attention in the 
literature. Their study was pioneered by Chernyak and Zhitnitsky (see |J for a comprehen- 
sive review). More recently, the results for p mesons were critically examined in fl2[]. In 
this paper we complete the update |12| by the reanalysis of SU(3) breaking corrections, see 
App. D 

A simple model of twist 2 distributions includes contributions of the three lowest con- 
formal spins ("S, P and D waves"): 



(u) 



Quu 



l + 3aU 



4 2^ 



(5.1) 



[u 



Quu 



l + 34£ + 4-(5£ 2 -i; 



(5.2) 



V 



We recall that u = 1 — u and £ = 2u — 1. For the K* meson distribution amplitude, 
designates the momentum fraction of a heavier (strange) quark. The values of the decay 
constants fy, fy and the Gegenbauer moments a%, a 2 are collected in Tabs. ^| and |] for 
V = p, K* and <p mesons. The scaling laws for the coefficients a\ and 4 are given in 
Eqs. ( [4.4q ) and (|3.58|) , respectively. 

The corresponding distributions, evaluated at p = 1 GeV, are shown in Fig. |I[ Note that 
the leading twist distributions for longitudinally and transversely polarized vector mesons 
are very similar to each other. 

In the following we neglect the masses of u and d quarks and do not account for p-oj 
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Figure 1: (a) Leading twist 2 distribution amplitude 4>\\(u) for a longitudinally polarized 
vector meson, (b) <j)±{u) for a transversely polarized one. Renormalization point is fi = 
1 GeV. 

mixing. In this approximation the couplings and distribution amplitudes of p^, p° and uo 
mesons are equal if one chooses properly normalized currents, i.e. (uu ± dd)/\/2 for the u 
and p° meson, respectively, in Eqs. ( PTFQ , ( g3g) , ( |2T2ip and ( gggj ). 



The model distribution amplitudes for the if* meson given by Chernyak and Zhitnitsky 
32, [3| involve an additional contribution ~ a 3 C^ 2 (^). We have not included this term since 



estimates of high partial waves from QCD sum rules are not reliable. Our estimates for aj- 
and differ significantly from the results of [32[ [§, see Ref. [12 for details. 



5.2 Three-particle distributions of twist three 

Three-particle twist 3 distribution amplitudes were defined in Sec. 2.3 and their conformal 
expansion is considered in detail in Sees. 3 and 4. Assuming SU(3) flavour symmetry, the 
lowest order contributions to the conformal expansion are 

V(a d , a u , a g ) = 5040 (a d - a u )a d a u a 2 , (5.3) 



9' 

1 



A{a d ,a u ,a g ) = 360a d a u a g \l +u 10 -(7a g - 3)], (5.4) 

T(a d , oc u , otg) = 5040 (a d - a u )a d a u a 2 g . (5.5) 

Our expressions for V and A agree with the corresponding "asymptotic distributions" in 
P, |33|. The result for T is new. An important point to note is that the contribution of 



leading conformal spin j = 7/2 to the distribution amplitudes V and T vanishes by virtue 
of G-parity invariance (in the SU(3) limit). Hence, if one takes into account the leading 
j = 7/2 contribution to the distribution A only, it is consistent to put V and T to zero; 
the expressions given in (|5.3|) and (|5.5|) correspond to contributions of the next-to-leading 
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conformal spin j = 9/2 and have to be taken into account together with the correction to 
A proportional to uj± , which has the same spin. 

The decay constants f^, f% and the few first coefficients u)£ k , ujj tk were estimated from 
QCD sum rules in Ref. 0, in particular 10 

/£ = (0.5 - 0.6) • 10- 2 GeV 2 , fl p = 0.2 • lO" 2 GeV 2 , = -2.1. (5.6) 

We expect that the accuracy for the last two entries is of order 30-50%. We have derived a 
new sum rule for fj p (see App. |CJ) from which we obtain the estimate 

fl p = (0.3 ± 0.3) ■ 10" 2 GeV 2 . (5.7) 

Again, the renormalization scale is fi = 1 GeV. The anomalous dimensions of the couplings 
can be found in Eqs. ( gjg§ ), ( P7Q and (|3~76| ). 



Estimates of quark mass corrections are difficult to obtain from such complicated sum 
rules. A detailed study of this point goes beyond the scope of this paper. For the present 
purpose we neglect SU(3) corrections to three-particle distributions and assume that the 
dimensionless couplings (,3 ,A ' T defined in ( |3.8| ) and ( J4.5| ) are the same for all vector mesons. 
The most interesting effect which we miss in this "poor-man's" approximation is that the 
leading conformal spin contribution ~ ?>Q0ctdCt u ct 2 g reappears in the distribution amplitudes 
V and T for K* mesons. These terms deserve a further study. Our preferred values for the 
parameters determining three-particle distributions are collected in Tab. [5| 



5.3 Two-particle distributions of twist three 



As repeatedly emphasized above, the equations of motion allow the elimination of two- 
particle distribution amplitudes of higher twist in favour of leading twist and three-particle 
distribution amplitudes as independent dynamical degress of freedom. With leading twist 
and three-particle distributions as specified above, we get the following exact expressions 
(including terms up to conformal spin 9/2): 



l + af^ (t^ + ^W -1 



Quu _ • j , ■ 6 
+ 3 5+ (3uu + u In u + u In u) + 3 <5_ (u In u — u In u) 



(5. 



,(*), 



u 



3£ 2 + ^(3£ 2 



v. 



4m 2 -3) + ^C 3 T (3-30£ 2 + 35£ 4 ) 



3 3 
+ -<5 + (l + £ \nu/u) + -5_£(2 + lnu + hiM) 



(5.9) 



10 The model distributions proposed in jlj include in addition the j = 11/2 terms with coefficients uj^ x = 
11.7, cj^q = 7 and u>2o = — 1-9- We do not include these contributions for simplicity and because the 
corresponding QCD sum rules are less reliable. 
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V 




K *± 


<\> 




0.032 


0.032 


0.032 


(J 


0.013 


0.013 


0.013 


S3 


0.024 


0.024 


0.024 




-2.1 


-2.1 


-2.1 


s + 





0.24 


0.46 







-0.24 










0.16 


0.33 







-0.16 






Table 5: Masses and couplings entering Eqs. (|5.8|)-( |5TTT|) . Renormalization point is fi 
1 GeV. We use m s (l GeV) = 150 MeV and put the u and d quark mass zero. 



g±\u) = Quu 



i + ak + {^4 + |c 3 A (i 



16 



OJ 



1,0 



35 
T 



C 3 y } (5£ 2 - 1] 



+ 6 5 + (3uu + u In u + u In u) + 6 <5_ (u In u — u In u) 



( 



9 



2 
105 



3 ~ 3 ~ 

+ - 5 + (2 + In u + In u) + - 5_ (2f + In u - In u) , 



(5.10) 



(5.11) 



where, for simplicity, we used asymptotic leading twist distribution amplitudes in the cor- 
rection terms proportional to quark masses ~ S±. u 

The resulting p meson distributions are plotted in Fig. ||| together with the corresponding 
asymptotic distributions and with the distributions calculated in the Wandzura-Wilczek 
approximation. It is seen that gluon corrections of twist 3 are generally important and tend 
to broaden the distributions. For the second moments we get (at the scale 1 GeV): 



[ l du (2u - l) 2 h^iu) = 0.24 (0.20), 
Jo 11 



du (2m - l) 2 hf{u) = 0.63 (0.60), 



(5.12) 
(5.13) 



n For realistic parameter values the correction in ai is small and can safely be neglected. 
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Figure 2: Two-particle twist 3 distribution amplitudes for the p meson. 

( l du{2u-lf g { l\u) = 0.25 (0.20), (5.14) 
Jo 

[ du {2u - l) 2 g { l\u) = 0.47 (0.40), (5.15) 
Jo 

where the numbers in parenthesis give the asymptotic values. As already mentioned, the 
oscillatory behavior of the distributions depicted in Fig. |2] is an artifact of the expansion 
in orthogonal polynomials and will be smoothened by contributions of higher-order partial 
waves. We expect, nevertheless, that our approximation is sufficient for calculating most 
overlap integrals that appear in physical applications. 

In Fig. |3| we compare the distribution amplitudes of p, K* and mesons, which differ 
due to the nonzero strange quark mass. Note that quark mass corrections to two-particle 
distributions in general involve logarithms of the momentum fraction and are not reduced 
to polynomials. This means that in this case the expansion in conformal partial waves does 
not correspond to an expansion in local operators, which is similar to what was observed in 
p3, for bilinear twist 4 operators. The quark mass effects are not large, but can result 
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0.2 0.4 0.6 0.8 1 
u 



0.2 0.4 0.6 0.8 1 
u 




Figure 3: Two-particle twist 3 distribution amplitudes for the p, K* and 4> meson. 

in a logarithmic enhancement of distributions close to the end-points u —>■ and u — > 1, 
see Eqs. fl5.8[)-( |5TTID and Fig. |3|. Because of that, the calculation of SU(3) breaking effects 
requires caution for physical observables which are sensitive to the end-point region where 
the linear approximation in m s breaks down.0 



6 Summary and Conclusions 

In the present paper, we have studied the twist three distribution amplitudes of vector 
mesons in QCD and expressed them in a model-independent way by a minimal number 



12 For finite quark masses, renormalization also gets complicated due to the absence of conformal symmetry. 
For example, Q k n _ k _ 2 could receive additional mixing not only with m q £}}} l _ 1 , but also with operators 



involving total derivatives (see Eqs. ( [3.67 ) and ( 4.47 )). Therefore, the scale dependence of /• 



fT. ,T 



-k-2 



may not be described by a simple extension of the anomalous dimension matrix involving the mixing 
with m q f v al l _ 1 , which is in contrast to the renormalization of twist 3 parton distributions. The complete 
clarification of this point is beyond the scope of this work. 
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of nonperturbative parameters. The one key ingredient in our approach was the use of 
the QCD equations of motion which allowed us to reveal the interrelations between the 
different distribution amplitudes of a given twist and to obtain exact integral representations 
for distribution amplitudes that are not dynamically independent. The other ingredient 
was the use of conformal expansion which, analogously to partial wave decomposition in 
quantum mechanics, allows one to separate transverse and longitudinal variables in the wave 
function: The dependence on transverse coordinates is represented as scale-dependence of 
the relevant operators and is governed by renormalization-group equations, the dependence 
on the longitudinal momentum fraction is described in terms of irreducible presentations of 
the corresponding symmetry group, the collinear conformal group SL(2,R). The conformal 
partial wave expansion is explicitly consistent with the equations of motion since the latter 
are not renormalized. The expansion thus makes maximum use of the symmetry of the 
theory in order to simplify the dynamics, which is related, in the perturbative domain, to 
renormalization properties of twist three operators. 

As it was known for some time [f2(], anomalous dimensions of twist three operators 
increase logarithmically with the spin. Like in the leading twist case, this property ensures 
convergence of the conformal expansion at sufficiently large scales and suggests that only 
the few lowest "harmonics" are important in calculations of physical observables. 

Based on this assumption, we have derived explicit and consistent models for all two- 
and three-particle distribution amplitudes of p, u, K* and mesons of twist two and three 
including contributions up to conformal spin j = 9/2. The relevant nonperturbative pa- 
rameters ( "reduced matrix elements" ) were estimated from QCD sum rules. The results are 
immediately applicable to a range of phenomenologically interesting processes like exclusive 
semileptonic or radiative B decays and hard electroproduction of vector mesons at HERA. 

Our formalism is — in principle — applicable to arbitrary twist, although its realization 
will become technically more involved. The application to twist four distribution amplitudes 
of vector mesons will be presented elsewhere. 
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A Formulae for Orthogonal Polynomials 

In this appendix, we collect useful formulae for the orthogonal polynomials which appear in 
the conformal expansion. 

Differentiation formula for Gegenbauer polynomials: 

4(1 - e)C 3 J 2 (0 = -{n + l)(n + 2)Cj/ 2 (fl. (A.l) 
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Differentiation formulae for Jacobi polynomials: 
d 



d 



(1 + 0^(0=2^ 



(n + 1) 



^(0 



(n + 2)(2n + 3) 



The equation ( |A.3| ) is obtained from ( |A.2| ) combined with (A.4) below. 
Recurrence formulae for Jacobi polynomials: 



ti+flpi 1 ' 1 '© 



(» + !)(» + 3) ^ + ^ + n+1 p ,„, K) 



(n + 2)(2n + 3)~ n ^ w ' " n w ' 2n + 3~ 

2(ra + 1) 
2n + 3 



(^ 0) (0 + ^?(0), 



^°' 0) (0 



n+1 
2n + 1 

n+1 
2n+ 1 



n -,(1,0), 



2n + 1 



TOO 



2n + 1 



pr\o+piV(o = (i+OPi ' 1 \o, 

Pi ' \0-P^(0 = (i-OPi m (0- 

Relations between Jacobi and Gegenbauer polynomials: 

(i + 0^ 0,1) (0 + (i - 0^ 1,0) (0 = 20^(0, 

(i + o^°' 1} (0 - (i - OPi lfi \o = acifico, 

(n + 2)pW)(0 = 2C^ 2 (0- 
Orthogonality relations for Appell polynomials ||19|| : 

(_!)*+» 



Jvaa d a u ap k ,(a d ,a u )J m , n (a d ,a u ) = 5 k+l , m+n 



2 k + l + 3 (k + I + 3)(2k + 21 + 5)\\ 



W, 



(A.2) 



(A.3) 



(A.4) 
(A.5) 
(A.6) 



(A.7) 

(A.8) 
(A.9) 

(A.10) 

(A.H) 
(A.12) 



(A.13) 
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where W, 



(k+i+i) 



d m+n Jk,i(ad, a u ) /da™da™ is a (A; + 1 + 1) x (A; + 1 + 1) symmetric matrix. 



This result can be obtained from the following relations: 



J Va a^ +l al +1 a 2 g J K i(a d} a u ) 







l) k+l k\l\ 



2 k+i+3(k + l + 3)(2k + 2l + 5)U 



(m + n < k + I) 
(m + n — k + I), 



(A.14) 



while the integral is in general nonzero for m + n > k + I. 
Integral formulae for Appell polynomials: 

d ru ru ]_ f \ 

— da d da u - \a d - h a u - 1 \ a d a u {l - a d - a u ) 2 J k ,i(o; d , a u ) 

du Jo Jo 1 — a d — a u \ oa d oa u J 



uu k\l\(-l) k ( k-l 



(1,1) 



2 (k + l + 2)\ \k + l + 3 



(A.15) 



— / da d / da u a d a u (l - a d - a u )J k ,i{a d , a u ) = 
du Jo Jo 

_uu k\l\(-l) k ( k-l N 

r k+l+2\S) r k+l+l\^>) 



2 (k + l + 3)\ \k + l + 3 



(A.16) 



d r u 

du 



[ da d [ da u - - \a d 

Jo Jo 1 — OL d — OL u \ 



_d_ 

da, 



- a,. 



a d a u {l - a d - a u ) 2 J k j(a d , a u ) 



uu k\l\(-l) k 
Y(k + l + 3)\ 



k + l + 3 



(A.17) 



r u , f u . 1 ( d d 

/ da d / da u - h - — 

Jo Jo l — a d — a u \oa d oa U/ 



a d a u (l - a d - a u ) 2 Jk,i{ot d , a u ) 



k\l\(-l) k 
(k + / + 3)!4 



-k + l 
2k + 21 + 5 



, k + l + 2 , (0i0) fn _p(o,o) fn 



(A.18) 



The results (|A.15|) — (|A.17|) can be obtained by differentiating and/or integrating the Appell 
polynomials Jk,i{<y d , ol u ) term by term. To obtain ( |A.18| ), it is convenient to calculate its 
derivative first, which can be done similarly to ( A.15|) , (|A.16| ) and ( A.17 ), and then integrate 
the result with the condition that it vanishes at u — 0. 
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B Conformal Expansion of Wandzura- Wilczek Contri- 
butions 



In this appendix we explain the structure of the conformal expansion of the Wandzura- 
Wilczek contributions to twist 3 chiral-odd two-particle distribution amplitudes, Eq. (|3.40| ). 
As mentioned in Sec. 3.2, the conformal spin assignment for these terms seemingly does not 
match the expansion in Eqs. ( |3.30|) and ( |3.31| ), which calls for an explanation. 



The basic idea is the following: the expansion derived in Eqs. (|3.30| ) and ( |3.31|) is based 
on the conformal expansion of the corresponding operators. The distribution amplitudes 
n» and are obtained as matrix elements of these operators between the vacuum and 
the longitudinally polarized p meson state. We call the state conformal and assign a cor- 
responding conformal spin, if it is annihilated by a conformal operator. The difficulty 
with the Wandzura- Wilczek terms is due to the fact that these contributions involve ma- 
trix elements over the p meson with a different (transverse) polarization. Working out the 
Wandzura- Wilczek contributions to fti , essentially corresponds to reexpressing these 
matrix elements in terms of similar matrix elements over the longitudinally polarized meson, 
using Lorentz symmetry. In our context it is important that the transversely polarized state 
is related (in the p meson rest frame) to the longitudinally polarized state by a spin rotation 
which does not commute with the generators of the collinear conformal group. Working out 
the necessary commutation relations, we reproduce the particular spin structure appearing 
in CT- 

It is convenient to work in the helicity basis: p meson states with A = ±1 correspond to 
transverse polarization, A = denotes the longitudinal polarization. Equation fl3.54|) reads 

yX«<o|n£|p(p,A = ±i)>, ( R1 ) 



we 



where the proportionality factor is irrelevant for what follows. Following Ohrndorf [le 
define a set of eigenstates \j, m) with conformal spin j and the "third projection" of the 
spin m = j, j + 1, j + 2, . . ., such that 

J 3 \j,m) =m\j } m), J 2 \j,m) = j(j - l)\j, m), 

J-\j,m) = -(j -m)\j,m- 1), J+\j, m) = (j + m)\j, m + 1). (B.2) 

The generators J±, J3 satisfy the canonical commutation relations of the algebra of the group 
of hyperbolic rotations, 0(2,1), and are related to the generators of the collinear conformal 
group by 

J + = -^=P, J_ = -j=K„ J 3 = l -(D + M*.), (B.3) 

where P M and are the usual generators of Poincare group, D is the generator of dilata- 
tions and generates conformal transformations. Note that J + and J_ are just "step-up" 
and "step-down" operators in this basis; the state with the lowest value of m, m m i n = j, is 
annihilated by J_. 
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Following the discussion in Ref. [jnj, it is easy to show that 

J_Q» = [J_, j£]|0) = 0, J s n£|0) = (n + 2)^|0>, J 2 ^|0) = (n + 2)(n + 1)Q», 

(B-4) 

so that we can identify 

|n + 2,n + 2) = ^|0), (B.5) 

and rewrite Eq. ( |B.1| ) as 

/Ja^ oc (n + 2, n + 2\p(P, A = ±1)). (B.6) 

This confirms that corresponds to conformal spin j = n + 2, as stated in the main text. 

To determine the corresponding contribution to hfj? , we have to recast Eq. ( P-6|) in 
a different form corresponding to a matrix element over the longitudinally polarized meson. 
For definiteness, take A = +1. In the p meson rest frame the A = +1 state is related to the 
A = state by the spin rotation 

\p(P = o, A = +1)) ex (M 23 + iM 31 )\p(P = 0, A = 0)), (B.7) 

where (M 23 + iM 31 ) is the step-up operator of ordinary angular momentum. \p(P, A)) is 
then obtained by a Lorentz boost in P 3 direction: 

\p{P,\))=U(u)\p(P = 0,\)), (B.8) 

where 

U{lo) = e-^ M ° S = e~ M -, th(w) = P 3 /P°. (B.9) 

We can thus write 

ff<£<x{*\p{P,\ = 0)), ( R1 °) 

where 

=WH(M 23 + zM 13 )W- 1 (w)|n + 2,n + 2). (B.ll) 

In the following we demonstrate that |^) is given by a superposition of three conformal 
states, 

I*) = C x \n + n + -) + C 2 |n + ~,n + -) + C 3 |n + ^,n+ -), (B.12) 

where the are C-numbers. If established, Eq. ( |B.12| ) shows that contributions of to 
the matrix element over a longitudinal p meson correspond to conformal spins j = n + § 
and j = n + |, which explains the pattern appearing in Eq. ( |3.40| ). 
To prove (B.12j ), we first note that 

U{uj) (M 23 + zMi 3 ) U' 1 ^) = M 23 ch(u) + M 20 sh(u) + i{M 13 ch(u) + M w sh(u;)}. (B.13) 

From this and Eq. flB.ll ), it follows 

= [exp(w)/2 (M 2 . + iMi.) - exp(-w) (M 2 * + zM 1# )] |n + 2, n + 2) (B.14) 
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with z» = (1,0,0,-1) and p^ = (1,0,0,1) (note that p ■ z = 2). Applying to (|KT4]) the 
commutation relations 

[J 3 , M 2 . + iMi.] = -- (M 2 . + zMl), [J 3 , ^2* + iMi,] = - (M 2 * + iM u ) (B.15) 

and 

[J-,M 2 . +iMi.] = 0, [J_,M 2 *+zMi»] oc K 2 + iK h [J_, K 2 + iHd] = 0, (B.16) 
it immediately follows that 

3 3 

(Ma. +iMi.)|n + 2,n + 2) oc |n+-,n + -), 

3 5 5 5 

(M 2 „ + zM u )|n + 2, n + 2) oc a x \n + -, n + -> + a 2 |n + -, n + -), (B.17) 

with C-number coefficients a^, which proves Eq. (|B.12| ). 

A similar discussion also explains the mismatch observed in Eq. ( |4.31| ) for chiral-even 
distribution amplitudes. 



C QCD Sum Rules for Expansion Coefficients of Dis- 
tribution Amplitudes 

The method of QCD sum rules in its application to distribution amplitudes of light mesons 
was pioneered by Chernyak and Zhitnitsky and is comprehensively discussed in |§. In this 
appendix we collect QCD sum rules and results for the twist 2 distribution amplitudes of the 
vector mesons p, K* and as well as for the twist 3 distribution amplitudes of the p meson. 
Numerical results presented below are obtained using the following input parameters: 

(150 ±50) MeV, (^G 2 ) = (0.012 ± 0.006) GeV 4 , 

(-240 ± 20) MeV 3 , (ss)(lGeV) = 0.8(gg)(l GeV), 

0.8(gg)(lGeV), (sagGs)(l GeV) = 0.8(qagGq)(l GeV), 

ct s (lGeV) = 0.56 (C.l) 

and assuming factorization of the vacuum expectation values of four-fermion operators. 
The SU(3) breaking effects in the sum rules are due to explicit corrections proportional 
to the quark masses, the difference in values of the condensates of strange and nonstrange 
quarks, and differences in the values of the continuum thresholds sq and Borel parameters 
M 2 . Instead of fitting the continuum thresholds separately for each meson and for each 
sum rule, in this paper we prefer to determine so, p from the simplest sum rules for vector 
(tensor) couplings and use the relations 

sq,k* - s ,p = m 2 K * - m 2 p , s o ,0 - s , p = mj - mj, (C.2) 



m s (lGeV) = 
(gg)(lGeV) = 
(qagGq)(lGeV) = 
A&L = 400 MeV 
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which are known to hold with reasonable accuracy for the spectra of resonances in the 
respective channels. Similar relations are assumed between the "working windows" in the 
Borel parameter. 



C.l Distribution amplitudes of twist two 

Most of the relevant formulas were previously obtained in [| |l2j. For the p meson, we 
quote the results from |T2|]. For the other mesons we present a new analysis which includes 
the radiative corrections calculated in [0 and the SU(3) breaking terms calculated in |32 . 
Unlike in Ref. R 1321], where sum rules for the moments of distribution amplitudes were 



derived, we prefer to consider the sum rules directly for the coefficients a n (Gegenbauer 
moments) in the conformal expansion, see for a discussion. 

QCD sum rules for even Gegenbauer moments can be derived from the diagonal corre- 
lation functions of the conformal operators introduced in Sec. 3.3 and Sec. 4.3: 



Di {±} =i /"Ae ira (0|raP>(y)O; il{±} (0)|0>, 



(C.3) 



see Eqs. ( |3.53| )-( p.55| ) and (|4.45[ )-( fl.47D for precise definitions 



One finds the following sum rules for vector and tensor couplings |30], |32| of the K* 
meson: 
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3M 4 



nr.i^fiGs) - ((qq) 2 + (ss)' 



81M 4 



In the sum rule for f3 — 9, ft\ — 64, and 7 



2tt 



12M 2 









Po 


v ^ 


Po) 


1 ot s 


G 2 ) + 


m s (ss) 


\ 71 


M 2 



X 



(C.5) 

310/9 (for three running flavours) is 



the two- loop anomalous dimension calculated in |p5 |. 

For arbitrary (even) Gegenbauer moments one obtains: 

3(n + l)(n + 2) f2 j\ riA „- m */M* 



2(2n + 3) 
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1 a 



2tt 2 



— M 2 (l — e~ s ° /M2 ) [ 1 duuuCl /2 (2u - 1) In 2 - + -L-m s (ss)(n + l)(n + 2) 
7r V /Jo u 2M Z 



I / a s 2 \ , . _ w I 



+ 24]^ \ ~ G j ( n + l)(n + 2) — m a (aagG 8 )n{n + l)(n + 2)(n + 3) 



2(2n + 3) 
1 a 



(qq)(ss)(n + l)(n + 2) + ^ ((gg> 2 + + l) 2 (n + 2) 2 , (C.6) 



g M 4 /V /V > glM 4 

3(n+l)(n + 2) T 



27T 2 7T 



M 2 (l - e ~ s o/ M2 ) £duuuCl /2 (2u - 1) (in u + Inu + In 2 | 



+ 24^2 ( V G ') + 3n - 2) - i -L_m.(^G?a)(n + l)(n + 2)(n 2 + 3n + 8) 

+ ((<7<7> 2 + (^) 2 )(^ - l)(n + l)(n + 2)(n + 4) + ^ m s (s S >(n + l)(n + 2). 

(C.7) 

The sum rules for p are obtained by setting m s zero and s to g in the condensates. For 
0, one has to replace q by s in the condensates and to double the terms linear in m s . On 
the left-hand sides one has to insert the proper meson masses. All renormalization scale 
dependent quantities are evaluated at the scale /i ~ 1 GeV. 

Gegenbauer moments with odd n are nonvanishing for the K* meson only. They can be 
determined most conveniently from the nondiagonal correlation functions 13 

NDp> = i |Ae w (0|raP>( 2/ )fii ±{ll} (0)|0). (C.8) 

The sum rules read [ 32fl : 

fK*(^)fK*m K *al(fi)e' m K' 

s > l2 ) + l(n + l)(n + 2)((ss)-(qq)) 
( n + if(n + 2) 2 {(sagGs) - (qagGq)), (C.9) 



3(n + l)(n + 2) f f iu]j . _ j _,„,.. ^ 
2(2n + 3) 

3 



, m a M i [l-e 

ilT 2 

1 



2AM 2 



3(n + l)(n + 2) f/ f()i/ . M1 , /l x, llW -,„i w-' 
2(2n + 3) 



13 Since pcrturbative contributions to D^ 71 ) vanish for odd n. Note also that from ND^ ^ one obtains the 

■T 
V ■ 



relative sign between fy and fT 
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m 



M 2 (l 



+ -(n+l)(n + 2)({ss)-(qq)) 



1 



48M 2 



(n + l)(n + 2)(n 2 + 3n + 4)((sac/Gs) - (qagGq)). 



(CIO) 



In these sum rules again the right-hand sides are to be evaluated at fixed scale p of order 
1 GeV. 

The results are collected in Tab. |], where the quoted errors are due to uncertainties 
in input parameters and to the variation of the Borel parameter within the range M 2 « 
(1-2) GeV 2 (for the p meson), with the value of the continuum threshold Sq„ = 1.5 GeV 2 
fitted to reproduce the experimental value of the vector coupling. As discussed in detail in 



the sum rule for the tensor couplings contains contaminating contributions of states 

which can be effectively taken into account by using a lower 
q = 1.2 GeV 2 . For other mesons, we assume validity 
of the relations (|G.2|). Note that in the numerical analysis of Gegenbauer moments we 



with the opposite parity 1 + , 
value of the continuum threshold s? 



substitute the couplings on the left-hand sides by their sum rules ( |C.4j ) and (|C.5| ) rather 
than using the values given in the table. 



C.2 Distribution amplitudes of twist three 

The vector /Z, and axial f£ twist 3 couplings are defined as the local matrix elements 



(0\d^ gG pX (iD.) - {iD.)gG 



u\p + ) 



gG pX (iD.) - {iD.)gG 



pX 



d\p-) 



v\^ P )fL + ■■■ 



(0\d~f^r,gG Xp u\p + ) = (0\uj^ 5 gG X pd\p ) = Pp(p P ei - p x e^ )ft + 



(C.ll) 
(C.12) 



and have been estimated using the sum rule approach in Ref. |3(| together with a few matrix 
elements of higher dimension (and conformal spin). The results are given in the textQ. 
The tensor coupling is defined as 



gG»(iD.)-(iD.)gG> 



u\p) = (ez)(pz) 2 m p fj p 



(C.13) 



and can be estimated from correlation functions of this operator with the vector and/or 
tensor current. The correlation function with the vector current is chirality- violating and is 
expanded in operators with odd dimension. The leading contribution of the mixed quark- 
gluon condensate, however, vanishes, and the first corrections comes from dimension 7 
operators whose vacuum expectation values are known only very poorly. For this reason 



14 Apparently the corresponding sum rules have never been published and are not available. We thank 
V. Chernyak for correspondence on this point. 
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this sum rule is essentially useless. From the "diagonal" correlation function with the tensor 
current, 



9 G^(iD.)-(iD.)gG> 



u(l/)u(0)<vd(0)|0), (C.14) 



d*ye^(0\Td(y)a, 

we obtain the sum rule 

% " /M2 m p f T U)flU) = [ S °dse- S / M2 + -(^G 2 ) + 7ra s (qq) 2 , (C.15) 

P Jp WJupW 7207r 3 J 36 \ 7T / 27M 2 m/ ' v ; 



which we have studied numerically. As a general feature of the sum rules for matrix elements 
of operators with high dimension, this sum rule is dominated at small M 2 ~(l-2) GeV 2 by 
the condensates of high dimension (four-quark operators, for the case at hand) and is not 
stable. At larger values of the Borel parameter the stability of the sum rule is very much 
improved, suggesting the value /^(l GeV) ~ 0.3 ■ lCT 2 GeV 2 . This number has to be 
considered as a rough estimate, however, since at large values of the Borel parameter the 
contributions of higher mass resonances and of the continuum are out of control. Note that, 
similar to the case of the twist 2 tensor coupling considered above, the sum rule (|C.15| ) 
includes contributions of states with opposite (positive) parity. Ascribing a 100% error to 
this result, we arrive at the range given in (|5.7p as our best estimate. 
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